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Introduction 


We have chosen to begin this course with a unit on recurrence relations 
because it is a topic we hope you will find fairly straightforward, while at the 
same time it is of both practical and theoretical value. Its practical value is 
partly that some mathematical models depend on recurrence relations and 
partly that numerical methods for solving differential equations depend on 
recurrence relations. The theoretical value of the topic is that it will 
introduce you to some mathematical ideas connected with the solution of 
linear equations, which will be followed up later in the course. 


In this unit recurrence relations will be examined in some detail, to see how 
they work, how to solve them, and how to cope with difficulties which can 
arise in using them. 


Study guide 


Sections 1, 3, 4 (television) and 6 (computing) should be studied in that order 
but Section 2 (tape) can be studied at any stage between Sections 1 and 5. 
This means that if you find yourself short of time before the TV programme 
you can leave your study of Section 2 until after the programme. To get the 
most out of the TV programme in Section 4 you should have studied 
Sections 1 and 3 beforehand. You can do the Exercises in Section 5 at 

any stage after studying Sections 1 to 4. (You may prefer to use these for 
revision later. ) 


Sections 3 and 4 both involve numerical work using a calculator. The results 
quoted in this unit were obtained using a calculator which displayed eight 
figures while working with ten figures. Your calculator will probably give » 
slightly different answers from those in the unit. 


It is planned that you should visit a terminal during your work for the first 
assignment as a means of practising some of the techniques and 
consolidating some of the ideas in Units 1 and 2. The package to be used for 
your computer work is described in Unit 2. When you make your visit to the 
terminal you should take both Units 1 and 2 with you. 


1 Introduction to recurrence relations 


If you have studied M101 or MS283 you will have met recurrence relations 
already under the heading ‘formula iteration’ and most of Subsection 1.1 may 
be familiar to you. If you think this is so, we suggest that you read it quickly 
and do the exercises at the end of the subsection. If you can do them, do not 
spend any more time on Subsection 1.1. 


1.1 What is a recurrence relation? 


A familiar question in an intelligence test is to be given a sequence of 
numbers and to be asked to determine the next number in the sequence. For 
example if you were given the sequence 1, 3, 9, 27, 81 you would probably be 
able to deduce that the next term (ie. member) in the sequence is three times 
the previous term and 3 x 81 = 243. The idea used is that the terms in the 
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sequence are related to each other. If we say that u, represents the rth term 
in the sequence and if we start with u,, the sequence up, u,, U5, U3,... above 
can be generated by the iteration formula 


Ur+4 aug 3u, (1) 


where we specify that u, = 1 is the initial term in the sequence. The formula 
(1) can be used to determine u, as 


u, = 3u, =3 x 1 =3. 

From this we use the formula (1) with r = 1, to give 
w= 3u, = 3 x3=9 

and so on, generating the sequence 1, 3, 9, 27, 81, 243,... 


The formula 
U, + ” ae 3u, 


is called a recurrence relation because it is a relation between succeeding 
members of the sequence which can be used recurrently to generate the 
sequence once the initial term is known. It is this idea that we want to 
explore in this unit as it is fundamental to many of the methods for solving 
problems on the calculator or the computer. 


The following examples illustrate some problems which give rise to 
recurrence relations. 


Example 1 

Consider the problem of determining how much money you will have in 
your deposit account after five years if you start with £1000 and the interest 
rate is 10% compounded annually. (Assume that you make no deposits or 
withdrawals in this period. ) 


After one year your £1000 will have produced an extra £100 in interest so 
that you now have £1100 in the account. In the second year you receive £110 
interest and the account increases to £1210. How do we express this problem 
mathematically? 


If we let u, represent the amount in the account after r years then we can 
write down u,,,, the amount at the end of the following year, as 


10 % 


interest 


amount at 
end of 


r th year 


amount at 
end of 


(r+i)st year 


= ° 
Ure { Ur = iu 

10° 
LC. Uy = iin. 


We start, at the beginning of the first year, with £1000 so that we write 
Ug = 1000. The calculator works very well on such problems, giving the 
following sequence. 


et tee ee 
1000 1100 1210 1331 1464.10 


After five years your account would have increased to £1610.51. 


1610.51 
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Example 2 

A more complicated banking problem would be one in which you start with 
£1000 in your deposit account and deposit £200 at the end of each 
subsequent year. If the interest remains at 10% throughout, the recurrence 
relation representing the growth of the account is 


amount at amount at 


end of end of deposit 
(r+ i)st year rth year 
U = Uv. + 10 Ue + 2.0 0 


u,, = Liu, + 200 


The following sequence was computed with the help of a calculator. 


tt 
1000 1300 1630 1993 2392.30 | 2831.53 


Your account, after five years, totals £2831.53. 


(In Section 4 recurrence relations will be applied to the similar but more 
realistic problem of mortgage repayments. ) 


So far we have only seen the case in which the recurrence relation gives u,, , 
in terms of u,. In the next example u,,, depends on the previous two terms. 


Example 3 

In 1202 an Italian merchant, Leonardo of Pisa, known to us as Fibonacci 
published a book on algebra, Liber Abaci (‘Book of the Abacus’), in which he 
posed the following problem: : 


‘How many pairs of rabbits will be produced in a year, beginning with a 
single pair, if in every month each pair bears a new pair which becomes 
productive from the second month on? 


The question is a little unclear as we do not know whether the first pair of 
rabbits is productive in the first month. We shall assume that they are not. 
We also assume that no rabbits die during the year. 


We let u, denote the number of pairs at the end of month r. Since there is 
one pair of rabbits at the beginning of the year we have uy = 1. Since we 
assume that this pair is non-productive there are no new rabbits in the first 
month and u, = 1. 


To set up the recurrence relation we need some way of determining the 
number of pairs at the end of month (r + 1). Since we have assumed that no 
rabbits die we can write the equation 


number of pairs _ | number of pairs number of new pairs born 
afterr+1months| | after r months in the (r + 1)st month 


ae of new pairs ae 


L€. u = u 
r+1 r in the (r + 1)st month 
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However, we know how many pairs were born in this month because each 
productive pair produces one new pair. The question tells us that the number 
of productive pairs are those pairs which are over one month old. Only the 
pairs of rabbits that were alive at the end of the (r — 1)st month satisfy this 
condition and there are u,_, of these. Hence the required recurrence relation 
is 

U,+4 = Uu, + Uu,_4- 
Using the known values u, = 1 and u, = 1 we compute the next terms as 

uo =u, +u,=14+1=2 

u,=U,+u, =24+1=3 


and so on. The following table gives the number of pairs of rabbits at the 
end of each month. 


No. of pairs. u, 


The solution to Fibonacci’s problem of determining how many rabbits there 
are at the end of the year is 233 pairs. 


The numbers in the above table are called the Fibonacci numbers. These 
numbers have some remarkable properties which have fascinated 
mathematicians over the years. In fact there is a current journal, the 
Fibonacci Quarterly, which is devoted to their study. We shall meet them 
again in Section 2. 


Example 4 3 

(Do not spend much time on this example as we shall not refer to it again.) 
If you have studied MS283 or M101 you will have met problems such as: 
determine the largest root of the cubic equation 


-=—20 4x + 1-6 


Dividing through by x” and taking all but the first term over to the right 
hand side of the equation gives 


(x + 1) 


es 


x=4- 


One method of attempting to locate the largest root is to use the recurrence 
relation : 


(x, + 1) 
x? 


r 


Xt =4- 


and to iterate in the hope that a point will be reached at which x,,, ~ x, 
(ie. X,4, 1 18 approximately equal to x,). The symbol ~ means ‘approxi- 


; = mately equal to’. 
If this happens an approximation to the root of the original equation has 


been found. 


Starting with x) = 3.5 the following sequence was obtained using the 


calculator. 


oo ee 
3.5 3.6326531 | 3.6489395 | 3.6508432 | 3.6510643 | 3.6510900 | 3.6510930 


\ 
| 
| 


3.6510934 


8 MST204 1 1.1 


Clearly this is converging to a limit, and we should be quite happy to say 
that the root is approximately x = 3.65109 to five decimal places. Why this 
sequence of numbers converges, why it converges to the largest root and to 
what accuracy the result can be quoted are questions which do not directly 
concern us here. The purpose of the example is simply to illustrate one of the 
ways in which recurrence relations can arise. 


In Examples 1, 2 and 4 each of the recurrence relations was of the form 


Un, =f (u,) (r =O 12.. J 
or yi = (n= Oz.) 


where f is a given function. It makes no difference whether we express the 
recurrence relation in terms of x, or u, since it is the relationship between the 
terms that we are interested in. Yet another way of writing the same 
recurrence relation is 


“uo =jiw4) ae ee ee 


We are interested in the relationship itself, not the way it is written down. 


The initial term that is needed to specify the sequence is called the initial 
condition. In Examples 1 and 2 we were given a value for u, and this initial 
condition enables us to compute as many terms in the sequence as we wish. 
There are other possibilities: you will frequently see, for example, an initial 
condition given for u,. This again is acceptable as we are still able to 
generate the sequence u,,U,,U3,U,... using the value for u,. 


The Fibonacci example is a little more complicated since each member of the 
sequence depends on the previous two and we needed two initial conditions 
to start the calculation. A general recurrence relation of this type is 


U,+4 = f(u,, ae 


where f(u,,u,_,) stands for some given formula containing one or both of 
the variables u, and u,_,. (In Example 3 we had f(u,,u,_,) =u, + u,_ 1.) 


It is also possible to have recurrence relations in which the formula for u,, , 
involves three or more of the preceding members of the sequence. This still 
satisfies the requirement that the terms in the sequence are related and that 
the formula can be used recurrently to generate the successive terms in the 
sequence. 


Finally, the formula giving the recurrence relation may depend on r as well 
as the previous terms. An example is the recurrence relation 
Uy + ov ru, 


Exercise 1 
Which of the following are recurrence relations? 


1 
(1) ae 5 (Ur + tr—1) 
(ii) u, = 3u, — 5 


a 1 4 
(iii) Xr+1 get xX, e2 x. 


2 a 
1 
(iv) Sa 
(v) u, = 3uU, +4 2e 2u, +2 


(Solution on p. 49). 
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Exercise 2 
In the following exercises, compute the next five terms. 


1 4 
(1) X45 == [s, - “| starting with x, = 1 
2 x, 
(ii) u.,, =Uu,+u__, starting with u, =1,u, =3 


(Solutions on p. 49). 


Exercise 3 

Find the recurrence relation which describes the situation in which £20,000 is 
deposited initially in a savings account that yields 10% interest per annum and £4000 
is withdrawn at the end of each year. When does the money run out? 


(Solution on p. 49). 


1.2. Classification of recurrence relations 


The use of recurrence relations is usually fairly straightforward. However, 
their analysis is more difficult, except for certain special types of recurrence 
relation and we are therefore going to restrict ourselves to these special cases. 
Before we can begin we need to introduce four new terms into our 
vocabulary in order to be able to distinguish the differences and similarities 
between recurrence relations. For example how would we compare the 
recurrence relation 


U, + 1 ak 3 U, + 2 
and the recurrence relation 


a 


X44 =4 
5 a 


r 


A recurrence relation of kth order is one where the difference 


between the highest and lowest subscripts is k. 


The order of a recurrence relation is important because the sequence 
generated by a kth order recurrence relation requires k initial conditions. 


In both the above examples the difference between the highest subscript 
(r + 1) and the lowest subscript (r) is (r + 1) — r = 1. Hence they are both 
first order recurrence relations. 


A linear recurrence relation is one that can be written as 


U,+4 = a,u, + DI 5 Cu, _ 2 = pagar, 2 P, 


where p, and the coefficients a,,b,,c,,... may depend on r, but 


do not depend on any of the ws. Otherwise it is said to be non- 
linear. 


The recurrence relations 
U4 ges 3 U, + 2 
and 
U, + 1 as ru, 
are linear while the recurrence relation 


x, +1 
Xo A x2 


r 


cannot be expressed in the above form and is thus non-linear. 
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A linear constant-coefficient recurrence relation is a linear 


recurrence relation where the coefficients a,,b,,c,... are 


j eee, 2h 


independent of r. (p, may depend on r.) 


This definition only applies to linear recurrence relations so we could not 
classify the formula 


in this way. 


The recurrence relations 
U4, = 3u, + 2 

and 
4o,, =3u, + 2P 

have constant coefficients while the recurrence relation 
Un+1 = MU, 


does not have constant coefficients because the coefficient of u, depends on r. 


A linear recurrence relation is homogeneous if p, = 0 for all 


values of r. Otherwise it is said to be nonhomogeneous. 


Again note that the definition only applies to linear recurrence relations. The 
recurrence relations 


Uy + — 3 U, 
and 
U,. + (=. ru, 
are homogeneous while the recurrence relations 
U,. + i 3 U, + 2 
and 
U., =r, +3 


are nonhomogeneous. 


Exercise 4 
Classify the following recurrence relations 


by completing the table: 


Constant- 
order? linear? coefficient? 
(i) u.,, = 1.1, + 2 
; (x, + 1) 
(ii) Xo, =4- So 


r 


Gi) 2 = 06 2 Oe. 


(iv) u.,, =4u,.,+7 

(v) “uo= 2Fu.., —u_; 

Vi) ye FF 
(vil) Uni 4 = U,U,_4 


(Solution on p. 50). 


homogeneous? 
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1.3. The solution of linear first order recurrence relations 


In certain simple cases it is possible to write down explicitly the general 

term in the sequence generated by a recurrence relation. The formula which 
does this is called the solution of the recurrence relation. In this subsection we 
shall examine how this can be done. The ability to write down solutions can 
be extremely useful in analysing the recurrence relation and it may also save 
some computation. 


The computation of the first few terms of the sequence generated by the 
recurrence relation 
U, + os 2u, 


starting with u, = 1 is fairly straightforward giving the following table of 
values. 


However, if we required the value of u.7, the computation would be 
extremely tedious if done, for example, using the calculator. In some cases it 
is possible to find a formula which gives the values of the sequence without 
reference to the recurrence relation. In this example you can verify that the 
formula 


u, = 2" 
yields the same values as the above table. If this is true for all values of n, we 
can determine u,,, as 


__ 4227 
Ux97 =2 


which is very easy to compute if you have a suitable calculator. This rather 
simple example can be generalized by considering the following questions: 


1. What would happen if we started with another value for uo, say Up = 5? 
If we start with uy = 5 in the recurrence relation 
U, + ee 2u, 


we obtain the following sequence. 


This is exactly the same as before except that each term in the sequence has 
been multiplied by 5. Thus we can write down our solution as 


u, =35 x 2" 


The effect of changing the initial condition from u, = 1 to uy = 5 is to 
multiply all the terms in the sequence by 5. 


2. Can we use these ideas to tackle the problem of solving the recurrence 
relation 
U, + i 2u, 


with uy = A, where A is any number? 
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If we start with u, = A and use the recurrence relation 
Met 2u, 


to generate the sequence, we obtain the following table. 


The pattern of the terms is emerging and we could confidently write down 
the general term 


“= AX 2 in = 01,2, 


By putting different values for A into this formula, we can arrive at the 
solution of the recurrence relation obtained earlier for the special initial 
conditions u, = 1 and u, = 5, and the solution for any other initial condition 
as well. A formula such as u, = A x 2” is called the general solution of the 
recurrence relation u,,, = 2u,; this distinguishes it from particular solutions 
such as u, = 2" or u, = 5 x 2", which can be obtained from the general 
solution by giving particular numerical values, such as 1 or 5, to A. The 
particular value of A is determined by the initial condition. 


Exercise 5 
Write down the first five terms in the sequence for the recurrence relation 


Uu,,4 = aU, 


r 


starting with the initial condition uy = A. Can you now write down the general 
solution for this recurrence relation? 
(Solution on p. 50) 


Exercise 6 

Write down the general solution of the recurrence relation 
Ur+4 ma Su,. 

What is the particular solution if u, = 3? 

(Solution on p. 50) 


Exercise 7 
(This is an important exercise which will be referred to in Section 3.) 
Consider the recurrence relation 


U,+4 es au, 


with u, = 1. For what values of a will the terms in the sequence 


(i) grow steadily 

(ii) diminish steadily 

(111) oscillate but grow larger 
(iv) oscillate and diminish 
(v) stay the same size? 


(Hint: Consider the following cases: a> 1,a=1,0<a<1l,a=0, -—1<a<JO, 
a= a< =f) 
(Solution on p. 50) 


14 The nonhomogeneous linear constant-coefficient first order 
recurrence relation 


We now turn to the slightly more difficult problem of finding the general 
solution of the recurrence relation 
Uy + es au, - Pp 


where a and p are given numbers which do not depend on r. This is the 
simplest kind of recurrence relation. 
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The most obvious way of trying to find the general solution is to start with 
Uy = A, calculate the next few terms using the recurrence relation and to 
look for a pattern. The first few terms are 
u;=auy + p=aA+p 
u, =au,+p=a(aA+p)+p=a’A+ap+p 
u; =au,+p=a(a’A+ap+p)+p=a@At+a*p+apt+p 
U4 = au; +p=a(ae°A+a*pt+apt+p)+p 
=a@A+a*p+a’p+ap+p 
After only the first four terms we can guess that the solution is going to be 
u,=a@"A+a" 'p+a" *p+---+ap+p 
=a"A+(a '+a"*7+---+a+4+1)p 
The expression in brackets is a geometric series. If we multiply it by (a — 1) 
we get 
Ba1@ +a? + et haeC4 e+ ae +a 
—q"-'—-+--—q*—a-1 
=q"— | 


(All terms cancel out except the first and last terms.) 


So that, provided a # 1, we have 


a" —1 
Set ee For a = 1 see Exercise 9. 


Hence we can express the solution as 


ia = 1) 


ae pie 


Finally, if we combine the terms in a”, we have 


For the general solution of the recurrence relation we regard A as an 
arbitrary constant. A particular value for A gives a particular solution. This 


last expression can be simplified further if we put B = A + 5 giving 


General solution of nonhomo- 
geneous linear constant-coefficient 
first order recurrence relation. 


where B is now regarded as the arbitrary constant which depends on the 
initial condition. Again a particular value for B gives a particular solution. 


The simplest particular solution of all is the one for B = 0. It is 


ee ee 
: ee 


if we choose the right initial condition, all the terms of the sequence will be 
the same. For general values of B, the term — p/(a — 1) still appears in the 
solution but there is now an additional term Ba" (see the formula in the box 
above). In Section 1.3 we saw that this term (with A instead of B) was 
precisely the general solution of the homogeneous recurrence relation 

Uu,+1 = au,. This new recurrence relation can be obtained from our original 


14 MST2041 1.4 


(nonhomogeneous) recurrence relation u,,, = au, + p by removing the 
nonhomogeneous term p (i.e. the term not containing any u). So the general 
solution shown in the previous box has the structure 


general 
solution of 


=| homogeneous 
problem 


boven | associated jawengy 


solution solution 


This structure is important; it will recur in other contexts. 


The structure we have just seen gives an alternative method for solving linear 
recurrence relations of the form 


U,+4 = au, +) 


which can also be used for finding solutions of higher order recurrence 
relations. We look for a particular solution and add it to the general solution 
of the associated homogeneous recurrence relation 


U4. 1 = au,. 


Provided a ¥ 1 this is easily done since there is a particular solution which is 
a constant. 


Example 5 
To find the general solution of the recurrence relation 


U,+4 = a 1 


we first find the general solution of the associated homogeneous recurrence 
relation 


Un+4 See 2u, 


obtained by dropping the nonhomogeneous term ‘—1’ from the original 
recurrence relation. This general solution is 


1=—f xz. 
We now look for a particular solution of 
i; = 26 4 


in the form u, = c where c is a constant. If such a solution exists it must 
satisfy the recurrence relation, Le. 


é= 2e =, 
so that c = 1. Therefore 


“oe i 


n 


is a particular solution. 


Using the structure of the solution, the general solution of the original 
recurrence relation is therefore 


oe ee = 


Exercise 8 
Find the general solutions of the following recurrence relations: 


(i) U, +4 sis Iu, 

(ii) a ., =118 2 2 
(ii) | ee 
(Solutions on p. 51) 


Structure of general solution 
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Exercise 9 Important Exercise 
What is the general solution of the recurrence relation 


Us 4 cee au, +p 


if a = 1? (Hint: Starting with uy = A look for a pattern in the sequence u,,u,, U3,... 
in using the recurrence relation 


U4. 1 = U, = p.) 
(Solution on p. 51) 


1.5 Finding particular solutions 


Once we have found the general solution of any recurrence relation, we can 
pick out the particular solution appropriate to any given initial condition by 
using this condition to evaluate the arbitrary constant, B. 


Example 6 
To find the particular solution with u, = 3 for the recurrence relation 
Uu,,, = 4u, + 6 


we substitute for a and p in the formula 


u, = Ba" — = t i 

to give 
u,=Bx 4-2. 

The particular solution can be found using the initial condition as 
uo = B-—2=3 

Le. 5. 


Hence the particular solution is 
ue=5x 4-2. 


If, instead of uy, we were given an initial condition such as u, = 30 in the 
above example we should find the correct value for B from 


mH =—Bx4 —2= 59 

ie. B = 2. The particular solution would now be 
u=2 x & —2. 

Exercise 10 


Find the particular solutions of the following recurrence relations with the given 
initial condition 


(i) lik, > = 148, tig = 1000 
(i1) 444 = 1.lu, + 200, iy = 1300 
(iii) as Sm, +B, ey = 5 


(Solution on p. 51) 


Exercise 11 

os ae 

a-—-1 ifaFl 
A+ np id= t 


Bd = 
Given that u, = 


express u, as a function of uo, a, n and p. 
(Solution on p. 51) 
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Summary of Section 1 


1. An equation relating a term in a sequence to one or more previous 
terms, which can be used to generate the whole sequence, is called a 
recurrence relation. 


2. A recurrence relation is of kth order if the difference between the highest 
and lowest subscript in it is k. The sequence generated by a kth order 
recurrence relation requires k initial conditions. 


3. A recurrence relation is linear if it can be written as 
U,+4 me aU, + ae = CU, — 2 a ee 2 P, 


where p, and the coefficients a,,b,,c,,... may depend on r, but do not depend 


gee) oii De 


on any of the w’s. Otherwise it is said to be non-linear. 


4. A linear recurrence relation is a linear constant-coefficient recurrence 
relation if none of the coefficients a,, b,, c, depend on r. (p, may depend on r.) 


5. A linear recurrence relation is homogeneous if p, = 0 for all values of r. 
Otherwise it is said to be nonhomogeneous. 


6. The general solution of a linear first order constant-coefficient recurrence 
relation 


uj. = au, +p 
Be = fae 
is aos yee 
A + np Ha=t 
where A and B are arbitrary constants. 


7. A particular solution for the recurrence relation is one in which the 
arbitrary constant is given a particular value. 


2 Linear second order recurrence 
relations (Tape Section) 


In the last section we gave an example of a linear second order recurrence 
relation when we discussed Fibonacci sequences in Example 3. The study of 
second order recurrence relations is probably more stimulating than first 
order relations because they give rise to more interesting sequences. The 
results of this section will also be useful when you study the units on second 
order differential equations and the numerical solution of ordinary 
differential equations. 


2.1 An example 
We begin with a modelling example taken from an economics textbook. 


Example 1 
(Do not spend too much time on the modelling process in this example.) 


We wish to examine the demands for the installation of telephone lines in a 
certain city. We define the following variables: 


Let p, = price for the installation of a telephone in year r 
q, = number of telephones installed in year r. 
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To produce a simple model we ignore, for example, the fact that there is a 
delay in installing the ’phone and we also ignore the effects of inflation. We 
also make the following assumptions: 


(i) If the price goes up there will be less demand for telephones. 


Suppose that, unknown to Telecom, the demand is given by a 
linear relationship of the form 


=e of 7 tees pe of ame 


installed in year r in year r 
Le. q, = 2000 — 3p. (1) 
(i1) If sales go up there will be more pressure on Telecom engineers who 


have to do the installation. To control this demand Telecom have to 
fix the price of installation, p,.,, at the beginning of year r+ 1. We 
suppose that the method used by Telecom to determine the increase 
in price is to make the increase equal to 10% of the expected 
increase in demand for telephones in that year. 


The price increase in year r + 1 is p,,, — P,. 


The latest figures for the demand at the beginning of year r + 1 will 
be qg, and qg,-;. Hence the latest estimate of the increase in demand 
per year is g, — g,_,. Taking the price increase to be 10% of this we 
have 

Pas P= ONG (2) 


From Equation (1) we can see that 


q, as q,—1 as ee = Pe) 


and hence we can substitute this into Equation (2) to obtain the linear 
second order recurrence relation 


Pr+1 es P, = 0.1(4, ex @2 3) 
= —O.3(p, oe P,-4) 


Le. p.., =0./p, + 0.3p__.,. 


In order to start the computation we need two initial conditions. Suppose 
that po, the price in year 0, is £100 and p,, the price in year 1, is £120; then 
we can determine the prices and sales in the following years using the 
recurrence relation. 


We have p, =0.7p, + 0.3po 
= 0.7 x 120+ 0.3 x 100 
= 114 
p, = 0.7p, + 0.3p, 
= 0.7 x 1144+ 0.3 x 120 
= 11538 


and so on. 
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The following table gives the prices and demands for telephones for the next 
few years. 


Year: No. of telephones: 
r g, = 2000: — 3p, 
0 1700 
1 1640 
z 1658 
3 1653 
4 1654 
2 1654 
6 1654 


(Prices have been quoted to the nearest penny and demands to the nearest 
integer. ) 


The numbers in the table appear to be settling down to a constant value. By 
the time you have finished this section you should be able to determine the 
general solution for this problem and should thus be able to confirm the 
behaviour of the sequence. 


One of the uses of this model could be to predict what might happen if, for 
some external reason, there was a sudden increase in price or if there were a 
sudden change, for example, in Telecom policy governing the method of 
determining the price increase in year r + 1. 


You can see that second order recurrence relations are only a little harder to 
use than first order recurrence relations. 


2.2 The solution of second order recurrence relations 


In the rest of this section we look at the solutions of a simple type of second 
order recurrence relation: the linear constant-coefficient homogeneous 
relation 


U,+4 a au, =< bu, _, 
where a and b are constants. 
By using the tape along with the following frames you should be able to 


discover for yourself how to obtain general solutions of some recurrence 
relations of this type. 


Now start the tape. 
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A \inear, second order, homogenéous recurrence relation 


Particular 
Solutions 


IS the general solution 
Un =A2™ + BIO™ P 


> 
8 
I 
\ 
> 
& 
a 
} 
Nh 
O 
> 
R 
ie 


Divide by 
Ax r-| 


a 
MW 
Ss 
8 

| 
» 
O 


A and B will not 
necessarily be 
Che same 


UW, =Ax2- or u =Bx10” 


Add solutions 


= nr 
A =O gives u, =Bxlo 


“ SAce” + Bx 10" 
B =Ogives u, =Ax2™ 
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Rule for solving linear, second order, homogeneous recurrence relations 
ee ee eee ee eee 


UW 


STEP ONE : Form the auxiliary equation 
Lo sax + b 


STEP TWO: Solve the auxiliary equation 


or xs = 44 


STEP THREE : Write down the general solution for 
the recurrence relation 


WwW ey 


nr 


Watch out ! 


There are snags still 
to come... 


Auxiliary equation with two equal roots 
ee ee 


STEP TWO : Solution of auxiliary equation: a 


Particular 
Solutions 
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General Solution of u_ | 


Solutions 
already found 


x ee = 
g Add Solutions ) 
Remember 
rm 16 not 


a constant! 
= (A + ie ae als 


Quxi liaury equati on 


2 different roots — 
2 equal roots -“ 
no roots? 


General rule 


STEP ONE: Auxiliary equation 


26 sees Oe eb 


STEP TWO: Solve auxiliary equation 


different roets equal 


A , Foots 


tro regl rects 


Solution is Solution 1s Solution depends 
ES n n zs n On Complex 
U, =ArX + Be Un =(A+ Bn DN mack ere. 
eee Unit 5 
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We can go slightly further than we did on the tape by looking at the 
auxiliary equation in more detail. When we solve the auxiliary equation 


x? —ax —b=0 
Do not Confuse this quadratic equation 
with the general quadratic 

ax +bx+C=20 
for which the formula method Qives 


J - Cc 


20a 


to find the two roots, using the formula method, we have 


_at./a? + 4b 


“e D 


and we have three cases to consider 


(i) a? + 4b> 0: the two solutions, 2 and y, of the quadratic will be different 
and the general solution for the recurrence relation is 


u, = Ad” + Bu" 


(ii) a* + 4b =0: the two solutions 1 and y will both be equal to a/2 and we 
have the general solution for the recurrence relation as 


u, = (A + Bn)2" 


(iii) a* + 4b <0: no real solutions exist. This case will be dealt with in the 
unit on complex numbers. 


Exercise 1 
Find the general solutions of the following linear homogeneous second order 
constant-coefficient recurrence relations 

(i) Un+4 — 7Uu, ee i2u =, (111) Ur+4 = U, + u,—4 

(ii) u,,=64—9u_, (NY) 44, —~ 03" Caz, 


(The third recurrence relation is the recurrence relation used to generate the 
Fibonacci Sequence and its solution will involve square roots.) 


(Solution on p. 52) 


2.3 Initial conditions 


For a second order recurrence relation the general solution involves two 
arbitrary constants A and B. The two initial conditions required to start the 
sequence can also be used to determine the values of A and B required for 
the particular solution of the recurrence relation. 


Example 2 
Find the particular solution of the recurrence relation 
“a =i A 
with u, = 5 and u, = 34. 
The general solution of this recurrence relation, which we found in the tape 
section, 1S 
“aAxKx PT +Bxei 


To determine the values of A and B we write down the equations for uy and 
u, in terms of A and B to obtain two equations. For this example we have 


“. = 2A + 108 34 (2) 


The solution to these equations can be found, for example, by subtracting 
twice the first equation from the second equation to give 


SB = 74 
Le. =. 


check in the 


Substitution in equation (1) now gives second equation | 
A+3=5 


1.€. A-co-Z. 


2x2 + 10x3 = 34+ 


MST20412.3 23 


Hence the particular solution is 
w= 2x2 + 3x10 


You may wonder whether it is always possible to find values for A and B 
given any initial conditions. Fortunately, it is possible to determine the 
arbitrary constants for any linear constant-coefficient homogeneous second 
order recurrence relation given the values of any two consecutive terms in 
the sequence. 


Exercise 2 
Find the particular solutions of the following recurrence relations satisfying the given 
initial conditions. 


re 64 = 74 — 12a, 2 aT 
ii) “#.,=64,— 9u_, 4 = 2,4, =7 
(iii) Uni, = U, Be U,— 4 Uy = 1, — 1 


liv) u.,, = 09n.-OC2a_, m= 28, =7 
(Solution on p. 52) 
Exercise 3 
For the general solution of the recurrence relation 
Ur44 as U, + U,—4 


given in Exercise 1(iii) what can you say about the behaviour of each part of the 
solution as n becomes large? 


(Solution on p. 52) 


Exercise 4 
For the recurrence relation 


Ra! 
U, +4 ae Qu, + Uu,_4- 


what happens to the term in the general solution as n becomes large? For the 
particular solution with up = 1, uw; = —% what happens as n becomes large? 


(Solution on p. 53) 


Exercise 5 
Find the general solution of the recurrence relation 


P,+1 = 0.7p, + 0.3p,_, 


used in modelling the telephones problem in Example 1. Find the particular solution 
if p) = 100 and p, = 120. What is the long-term behaviour of this particular solution? 


(Solution on p. 53) 


Summary of Section 2 


1. Provided two initial conditions are given we can generate a sequence of 
numbers using a second order recurrence relation. 


2. To obtain the general solution of the linear homogeneous constant- 
coefficient second order recurrence relation 


U.+4 = au, - bu, 
we first consider the auxiliary equation 
ae 
x“ =ax+b. 
There are 3 cases to consider 


(i) a*+4b> 0: the auxiliary equation has two distinct solutions, 
A and p, and the general solution for the 
recurrence relation is 


u, = Ad" + By 


Important Exercise 
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(ii) a*+4b=0: the auxiliary equation has two equal solutions 
and the general solution is 


u, = (A + Bn)A" 


(iii) a? + 4b <0: the auxiliary equation has no real solutions. 
(See Unit 5.) 


3. To obtain a particular solution of the linear homogeneous 
constant-coefficient second order recurrence relation we use the two 
initial conditions to evaluate the arbitrary constants A and B in the 
general solution. 


3 Numerical difficulties in using 
recurrence relations 


3.1. Errors in computation and data 


When the solution to a problem involves a lot of calculation, using either a 
calculator or a computer, we have to be aware that we are not going to get 
the exact solution to the problem. There are three distinct types of error that 
may affect the solution: 


(i) Rounding errors, due to the fact that numbers cannot be stored exactly, 
may introduce a significant error in the solution. 


(ii) The data given may not be exact. 


(iii) The mathematical formulation of a problem may involve 
approximations or simplifications. 


In this section of the unit we shall only be looking at the problems which 
arise due to the first two types of error. We begin our discussion by 
describing two examples which illustrate some of the difficulties in using 
recurrence relations. The first is a very simple first order recurrence relation 
with constant coefficients while the second example is more realistic. By 
looking at the difficulties encountered using constant-coefficient recurrence 
relations you should see why the second problem causes trouble. 


Example 1 
Use the recurrence relation 
“u,, = 10a, —3 


starting with uo = 4 to evaluate u,, on your calculator. Compare your answer 
with the true solution. 


The only difficulty I have is in storing } in my calculator which displays 8 
figures but actually calculates using 10 figures. Here are the displayed results 
I obtained. 


0 0.333 33335 

1 0.333 333 33 0.333 
2 (1333 35555 0.33 
3 0.333 333 3 0.3 

4 0.333 333 0. 

5 0.333 33 


You may get slightly different results if you use a different calculator. 
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The general solution for the recurrence relation, using the formula developed in 
Section 1, is 


u,=Bx 10" +, 


and the particular solution required is u, = 4. My computed solution is 
U;2 = — 33 whilst the theoretical solution is u,;, = 4. Hence my computed 
solution is totally wrong. 


Example 2 
The values of the integral 


1 
ee | Me dc 
0) 

are required for various values of n. One way to do this is by numerical 
integration, given a particular value for n. However, integrating by parts once 
gives 

1 

i= ($e | oe a 
0) 
1 
= =n| ae ee 


0 


1 

But | x"~* e*~! dx is I,_, so that we obtain the recurrence relation (using r 
0) 

instead of n) 


Pei -7i, 


What we have done is to replace the explicit definition of J,, in a form we 
found difficult to handle, by a first order recurrence relation, which we can 
use to generate all the values for the integrals. Note that it is not a constant- 
coefficient recurrence relation. 


To start the sequence we need an initial condition. We note that I, is given 
by 


1 
= | gt dx = [e"*]}}=1—e ! =0.632 120558828 57... 
O 
I used my calculator to compute the first 14 terms starting with 


Ig = 1 — exp(—1) = 0.63212056. The numbers produced appear in the 
following table. 


) 0.63212056 0.10094291 
1 0.36787944 0.09151379 
2 0.26424112 0.08486208 
3 0.20727665 0.06651712 
a 0.17089342 0.20179456 
2 0.14553291 —~ 1, 0235293 
6 0.12680255 23.72661 

f 0.11238214 


You may get slightly different results if you use a different calculator. 


Now should we accept these results? Is the value of I,, approximately 23.7? 
To see that I have produced silly answers for I,, and I,, we return to the 
original definition of I, as an integral. Since the value of I, represents the 
area under the graph of y, = x"e*~* over the interval [0,1], it is useful to 
look at some of these graphs. 
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0 02 04 06 16 30x 0 O02. 04.06 US 29-2 


0.5 


0 02 04 06-02 10.. x 0 U2 Us 06 028 10. -& 


Figure 1 


In each of the graphs you can see that the function we are integrating is 
positive in the interval [0,1] and we know that ah hg) eg 
point in the interval [0,1]. We conclude from this that the integral must be 
positive and also the values of the integrals should decrease as n increases. 
Our knowledge of the problem indicates to us that the values I obtained for 
I,5, 1,3; and I,, cannot be anywhere near the correct values. 


This is a very interesting problem because we know how the solution should 
behave. For many practical problems we will have this insight and it can 
often tell us that things have gone wrong. The sole cause of the trouble in 
both the above examples is an apparently insignificant error in the initial 
condition. In both examples I made no errors at all in using the recurrence 
relation. 


Both of the above examples illustrate a phenomenon called ill-conditioning. A 
problem is said to be ill-conditioned if a small change in the data for the 
problem results in a significant change in the solution. This definition is a bit 
vague because we need to know what the ‘data’ are and also what we mean 
by ‘changes in the data’ and ‘significant changes in the solution’. We shall 
illustrate what this means by looking at the problem of finding uy for a large 
value of N using the linear first order recurrence relation 


U,+4 ee au, + P 


with u, given and a and p known constants. 


The data for the above problem are the values of u,, a and p. Each of them 
may be the result of some measurement or modelling process. As such they 
may be subject to measurement or modelling errors so that even before we 
start our computation we do not have exact information. 


From Exercise 11 in Section 1 we can write down the particular solution for 
the above recurrence relation as 


wt Polar Ee fal 


n 


Uy + np ita = 1 
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This formula tells us that the solution u, depends on each of the three data 
variables uy, a and p so that a change in any of these variables will result in 
a change in the solution. To do a complete analysis we would have to 
determine how u, varies for small changes in each of the three variables but, 
for simplicity, we are going to concentrate on one special case: ie. How do 
small changes in uy affect u,? 


Exercise 1 
Find the solution, u,, of 


u.,, = 10u, — 3 


given that u) = 3 + ¢, where ¢ is a given (small) number. 


Does a small change in uy lead to a smaller or larger change in iat 


(Solution on p. 53) 


3.2 The effect of absolute errors in the initial condition 


Examples 1 and 2 in the last subsection were both examples of small errors 
in the initial condition u, causing catastrophic errors in the solution. Here 
we assume that the coefficients in the recurrence relation are given exactly 
and that no rounding errors are made in any of the computations. The only 
error introduced will be an error in the initial condition. 


Before we analyse the propagation of errors we shall briefly review the ways 
in which errors are measured. The error in an approximation xX to a quantity 
whose exact value is x can be expressed in two ways: 


(1) The absolute error in the approximation is defined as x — x 
ie 


(ii) The relative error in the approximation is defined as 


It depends on the circumstances which is the appropriate measure of the 
error to use. In this subsection we examine the propagation of absolute 
errors while relative errors will be discussed in more detail in the next 
subsection. 


Example 3 

The value of z is to be approximated by 3.142. Since the correct value of z is 
3.14159265... we can estimate the relative and absolute errors in the 
approximation. | 


The absolute error: X — x = 3.142 — 3.14159265 ... 
= 0.00040734... 
~ 0.4073 x 107? 


X—x  0,00040734... 
~ 3.14159265... 
= (,0001296 ... 


= 0.1296 x 107° 


The relative error: 


To illustrate the way errors can build up in a recurrence relation calculation 
we go back to Example 1, which was to calculate u,, given that u, = + and 
u,+, = 10u, — 3. Let us write uy, u,,u5,... for the solutions of the recurrence 
relation, starting from uy = $ and %),U,,U5,... for the solution starting from 
a slightly different value, 4) = 4 + «. By Exercise 1, the solution starting from 
3+ is \ 


1 
to te x 10! 
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How does the error in %, depend on n? With the initial condition uy = + the 
nth term in the sequence is u, = 3. Hence the error in u, is given by 


1 1 
i= 8, = |= +e x 10] —= 
,-4= G+ex 10) 3 
=e 
= (iy — Uy) x 10° 
ie. the absolute error in H, is 10” times the error in %. We call the number 
10" the scale factor for the computation. We can see now what went wrong Scale factors are treated in M101 
with the sequence generated in Example 1. By the time we compute u,, the and MS283. 


error in i, has been magnified 10' times and even with a calculator which 
works to 10 figures any small absolute error in uy will swamp the solution. 


To describe such situations we introduce the term absolutely ill-conditioned, 
defined as follows. 


A problem is said to be absolutely ill-conditioned if a small 
absolute error in the data gives rise to a significantly larger 


absolute error in the solution. If the absolute error in the 
solution is smaller than the absolute error in the data the 
problem is said to be absolutely well-conditioned. 


In the example we have just considered, a small error ¢ in the data produces 
an error 10/2 times as large in the solution, u,,, so this problem is absolutely 
ill-conditioned. 


Using this definition of absolute ill-conditioning we shall now look at the 
general linear constant-coefficient first order recurrence relation 
Uy + es au, - Pp 


where the problem is to determine u, given an initial value uo. If the 
approximation to uy is Uy = uy + & the approximations U,,U,,... to the rest 
of the sequence will (assuming all calculations are done exactly) satisfy 


U,+4 = au, + Pp 


with the initial condition Up = Uo + €. 


Subtracting the two recurrence relations we have 


Su. is to be 


treated as a single 
Symbol. It does 


Not mean § x UW. 


Uy +1 Se U,+4 = a(u, = Uu,). 


Defining du,=u, — u, as the absolute error in u, we obtain the 
homogeneous recurrence relation for du, as 


Ott, 5 = aou, 
We have already examined the behaviour of this recurrence relation in 
Exercise 7 of Section 1 and we can summarize those results as 


(i) Ifa>1ora< —1, du, will grow in magnitude as n increases 
Gi) If —1<a<1, du, will not grow in magnitude as n increases. 


We can also work out the scale factor for the computation of u, since the 
solution of 

Ou, + —— aou, 
iS Ou, =U, — u, = a" Ou 
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Hence the scale factor is a". 


For large values of n we have the result that: 


The problem of using the recurrence relation 


Up», = au, + P 


to determine u,, for large values of n, given up, is absolutely ill- 
conditioned with respect to small changes in uy if |a| > 1. If 
|a| < 1 the problem is well-conditioned. 


Exercise 2 
Determine the scale factors and comment on the absolute conditioning of the 
following problems 


(i) Find u,, 
using u.,, = 3u,+2 
with uy = 1 

(ii) Find u,, 
using u,,, = —3u, + 6 
with u, = 1 


(Solution on p. 53) 
We can extend our analysis to Example 2 where we are using the recurrence 
relation 

i, a 1 ae ae 1 


to compute the values of the integrals using J, = 1 — e~*. If we make a small 
absolute error in J, then the recurrence relation is 


I.=1-—rl2, 
Subtracting the two equations gives 
I, —1, = —r(l,_, —I,_-;) 
and putting 6], = I, — I, to denote the error in I, we have 
ol, = —r ol,_, 


Hence any absolute error in I,_, will be magnified by —r to give the 
error in J,. To compute the scale factor for I,, we see that 


of = —nol_; 


n 


= n(n — 1)oI,_, (Since 6f,_, = —(n — 1)0I,_,) 


= (—1)'n! 51, 


Hence the scale factor for this problem is (—1)’n! This becomes large very 
rapidly as n increases. For example if n = 12 then the scale factor is 12! 

= 479 001 600. We can conclude that the problem in Example 2 is absolutely 
ill-conditioned. 


Generally, the problem of determining u, given u, using the recurrence 
relation 

U,+14 aay aU, + P, 
will be absolutely ill-conditioned with respect to small changes in u, if 


la,| > 1 for all values of r. 


Exercise 3 
We wish to use the recurrence relation 


U,+ 4 = ru, — 1 
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with up = 1, to compute u, for large n. Determine the scale factor for the problem and 
comment on the absolute conditioning of the problem. Does the absolute 
conditioning depend on the initial condition up? 


(Solution on p. 53) 


Exercise 4 
Use the recurrence relation 
u., = Ka, — 35 
to compute u,, for each of the three initial conditions 
(i) uy =1 
(ii) , = 0.9999 
(ii1) Uy = 1.0001 
Comment on the answers obtained. 


(Solution on p. 53) 


3.3. The effect of relative errors in the initial condition 


The results of the last exercise should have convinced you that absolute 
errors are not necessarily the best measure of the errors in the computation. 
If the terms in the sequence are also growing rapidly the growth of the 
absolute error may not have a significant effect on the value of u,. Changing 
uy from 1 to 1.0001 changes the exact value of u,, from 6 666 666 667 to 

6 667 666 667 (i.e. the solution has changed by one million). However the 
value of uo is so large that the magnitude of the absolute error might be 
considered acceptable. 


Look at what happens to the relative errors in this example. Writing uy) = 1 
and Ui, = 1.0001 we have u,,. = 6 666 666 667 and wu, = 6 667 666 66/7. 


a 0.0001 
The relative error in Up, is = Z a ee 7 x 10 
0) 
hg = 1 000 000 
The relative error in U9 1s —— = 6666666667 > 0.00015 
= 0.15 x 10°° 


The magnitude of the relative error has grown by a factor of only 1.5 and, in 
situations where relative errors are more important than absolute errors, this 
might be deemed acceptable. 


To emphasize the distinction between the growth of relative and absolute 
errors we define the term relatively ill-conditioned as follows: 


A problem is said to be relatively ill-conditioned if a small 
relative error in the data gives rise to a significantly larger 


relative error in the solution. If the relative error in the 
solution is smaller than the relative error in the data the 
problem is said to be relatively well-conditioned. 


Relative ill-conditioning occurs far less frequently than absolute ill- 
conditioning. Examples 1 and 2 in Subsection 3.1 were, however, illustrations 
of relative ill-conditioning. The feature which was common to both examples 
was that the scale factor was large while the solution we required was small. 
For first order recurrence relations these are precisely the conditions required 
for relative ill-conditioning. 
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To see how relative ill-conditioning can occur, we examine the problem of 
using the first order recurrence relation 


U,41 = au, + P 


to determine u, given up. 


We have seen in Section 3.1 that the general solution can be expressed as 


a—1l 


Uy + Jo ifaA#Al 
Uy + np fa=1 


Now we know that absolute errors grow if |a| > 1. Suppose that we do have 
|a| > 1. Then u, is given by 


a—l 


The first term on the right hand side is increasing rapidly in magnitude and 
normally the size of u, will increase rapidly as n increases. However for one 
initial condition the value of u, will remain small, i.e. when 


Pp 
a—1 


‘i= 


in which case the solution is 


For the special initial condition uy = -—, therefore, the true solution 
a eS 


remains small while the computed solution grows rapidly: the problem is 
relatively ill-conditioned. It turns out that relative ill-conditioning also occurs 
when ug is close to the special value — = without being exactly equal to 


it. The analysis of this case will not be given here but you should know that 
it occurs. 


These results can be summarized as follows: 


The problem of determining u, for sufficiently large n using the 
recurrence relation 


U,+4 = au, Pp 


is relatively ill-conditioned with respect to small changes in uy 


if 
(i) lal >1 


P 


and a, = So 


Exercise 5 
For the recurrence relation 


U,+4 es = su, + 2 


show that the problem of determining u,, given u, = 4 is relatively ill-conditioned, 
with respect to small changes in up. 


(Solution on p. 53) 
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Exercise 6 
Find the values of u,,. and U,, for the recurrence relation 
1 
es aes ar +3 


with the initial conditions 

(i) Uy = 

(11) Uy =6+6 

What are the absolute and relative errors in U,,? Compare these with the absolute 
and relative errors in uy). Do your findings agree with the results on absolute and 


relative ill-conditioning? 


(Solution on p. 54) 


Exercise 7 
(i) Use your calculator to compute the next five terms in the sequence using the 
recurrence relation 


U, + 1 ms 50u, Gx = 
if the initial condition u, = 4 is to be approximated by u, = 0.142857. 


(ii) Compute the next five terms using the same recurrence relation if the initial 
condition is u, = 4, approximated by u,) = 0.166667. Comment on the results of (i) 
and (ii). 


(Solution on p. 54) 


We can extend the above argument to non-constant-coefficient linear first 
order recurrence relations similar to the one satisfied by the integrals in 
Example 2. In that example, as in the constant-coefficient case, there is 
precisely one initial condition for which I,, remains small as n increases. Any 
other initial condition from I, = 1 — e~* will result in values of I, which 
increase rapidly as n increases. 


Generally, the problem of determining u, for sufficiently large n, given the 
recurrence relation 

Ur+4 es a,u, = Pr 
and the value of uy, will be relatively ill-conditioned, with respect to small 
changes in up, if 
(1) |a,| > 1 for all values of r 


(i1) u, is close to an initial condition for which the solution remains 
small (i.e. grows less rapidly than the solution of the associated 
homogeneous problem). 


Exercise 8 
In this question we look at the recurrence relation 


Uy, =1+ (r+ 1)u, 
Compute the first 11 terms with 
(i) Uy = —1.7182818 
(ii) uy = —1.7182819. 
What conclusions can you draw from these results? 


(Solution on p. 54) 


Summary of Section 3 


1. Errors arise in numerical calculations which may or may not severely 
affect the solution. The absolute error in the computed solution u, compared 
with the true solution u, is u, — u,. The relative error in the computed 
solution u, compared with the true solution u, is given by (U, — U,)/Up. 
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2. A problem is said to be absolutely ill-conditioned if a small absolute error 
in the data gives rise to a significantly larger absolute error in the solution. 
A problem is said to be relatively ill-conditioned if a small relative error in 
the data gives rise to a significantly larger relative error in the solution. The 
growth of relative errors is particularly relevant when dealing with 
computations on the computer or calculator. 


3. In this section we have concentrated on one of the more obvious cases in 
which relative and absolute ill-conditioning arises. Whilst relative ill- 
conditioning is a comparatively rare phenomenon you should be aware of 
the special cases in which it arises as described in this section. 


4 For the recurrence relation 
U, + 1 — au, + Pp 


with a given initial condition u, the scale factor for the computation of u, is 
a". If |a"| is significantly larger than 1 the problem of computing u,, is 
absolutely ill-conditioned with respect to small changes in uy. For sufficiently 
large n this occurs whenever |a| > 1. Similar results hold for the problem of 
computing u, using the more general recurrence relation 


Up+4 a8 aU, = P, 
which is absolutely ill-conditioned with respect to small changes in up, for 
sufficiently large n, if |a,| > 1 for all values of r. 
5. The computation of u, using the recurrence relation 

U, + , wees au, “3 p 


is relatively ill-conditioned with respect to small changes in u, for sufficiently 
large nif 


(i) lal >1 


(11) Uy xX — 
The problem of computing u, using the more general recurrence relation 


Un+4 = a,u, ar P, 


is relatively ill-conditioned with respect to small changes in uo, for sufficiently 
large n, if 


(1) |a,| > 1 for all r 
(ii) Uy is close to an initial condition for which the solution remains 
small. 


6. The following flowchart summarizes our findings for constant-coefficient 
linear first order recurrence relations. 
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PROBLEM 
Find Un for large 
values of n using 
Mr+j = 24+ p 
With uw, given. 


Problem is both 
relatively and absolutely 
well- conditioned with 

respect to small changes 


in We 


TRUE 


FALSE 


Problem is absolutely 
ill- conditioned with 

respect to small changes 
In Uw, 


- Problem is relatively 
well - conditioned with 
respect to small changes 
In We 


Problem is relatively 
ilk- conditioned with 
respeck to small 
changes in U, 


FALSE TRUE 


4 Living with ill-conditioning 


This section consists of two parts. The first is based on a television 
programme about a specific problem where ill-conditioning arises in a 
modelling context; the second (not based on the television programme) shows 
how to cope with ill-conditioning under certain circumstances. 


4.1 Modelling a mortgage (Television Subsection) 


Read these notes before viewing the programme. 


The story of the programme is a fairly straightforward analysis of the TV1 
problem of determining the amount to repay each month after taking out a 

mortgage. It is not difficult to express this problem in a form which leads 

naturally to a recurrence relation. 


However the purpose of the programme is not to teach you how to solve 
recurrence relations of this type. (This is done in Section 1 and you are asked 
to solve the recurrence relation, derived in the programme, as an exercise to 
be done before viewing the programme.) Rather the programme looks at two 
topics involved in the solution and these two topics form the two halves of 
the programme. 


1. What are the modelling steps to be carried out in formulating the 
problem as a recurrence relation? 


2. How does ill-conditioning show up in the solution of the mortgage 
problem? 
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Before watching the programme you should do the following exercise. The 
solutions can be found in the programme notes which follow. 


Exercise 
(i) Find the general solution of the recurrence relation 


u,,=Lita,-UsM 
where M is a constant. 
(ii) Given that M = 258, for what values of u, is the problem 


(a) absolutely ill-conditioned? 
(b) relatively ill-conditioned? 


Now watch the television programme ‘Modelling a mortgage’. 
Read these notes after viewing the programme. 


One of the most important financial decisions in many people’s lives is the 
decision whether or not to take out a mortgage to buy a house. In this 
programme we consider the case of a £20000 loan which is to be paid back 
over 25 years when the interest rate is 15° per annum. 


In accord with one of the basic principles of modelling the problem is kept 
as simple as possible, at least in the first stages of the modelling process. No 
account is taken of any extra payments, such as house insurance, which can 
be added on to the monthly repayment. We simply have a mortgage which is 
to be gradually repaid over a period of 25 years. 


The way in which the building societies do their calculations is as follows. At 
the beginning of each year the 15 % interest is added to the amount 
outstanding. During the following 12 months the repayments are made at 
monthly intervals. For example, the interest charged in the first year is 15% 
of £20000 which is £3000. If M is the monthly repayment, the amount owed 
at the end of the first year is 


20000 + 3000 — 12M 


The building society will then charge interest in the second year, based on 
this amount. As this happens each year we can write down our basic 
equation as 
Amount owed 
next year 


pee owed 


+ [interest] — [repayments | 
this year 


Suppose u, represents the amount owed after r years. Our equation can be 
written as a recurrence relation: 


15 
ue 10M 
1 U, + 100" 


U, + 
Le. u,,, =1.15u,— 12M (1) 


This linear, first order, nonhomogeneous recurrence relation is the one you 
solved before watching the programme. Its general solution (see the box on 
p 13) 8 


u, = B(1.15)” + 80M (2) 
where B is an arbitrary constant. 
This solution contains two unknowns B and M and we need two conditions 
to evaluate them. The two extra pieces of information which we have are: 


(i) The initial mortgage is £20000, Le. uy = 20000. 
(ii) The mortgage is repaid in 25 years, Le. u,, = 0. 
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Putting these values into the general solution (2) gives two simultaneous 
equations: 
20000 = B + 80M 
0 = B(1.15)?? + 80M 
Subtracting the second equation from the first we have 
20000 = B(1 — (1.15)?*) 
and so 
B = —626.58696... 
Substituting this value into the first equation gives 
20000 — B 
M = oe ee: 
= 257.8324... 


The monthly repayment figure M is not an amount we can pay exactly. The 
best we could do is to pay £257.83 per month or in whole pounds, as in the 
programme, £258. Using M = 258 we obtain the recurrence relation 


uu, , = 1.15u, — 3096 


This recurrence relation could now be used to determine how much you 
would owe each year. Rounding off each value of u,, as it is calculated, to the 
nearest pound we obtain the following table of values: 


CONNMN BR WN KF © 


A rounding error of 17 pence per month plus the rounding of the interest 
payments to the nearest pound has resulted in the payment of £427 too 
much. Now 12 repayments of 17p is £2.04. Thus over 25 years you might 
have expected this to affect the result only by about £51. However, the pre- 
programme exercise should have prepared you for the magnification of this 
error. This problem is absolutely ill-conditioned for any value of uy since the 
factor 1.15 multiplying u, in the recurrence relation is greater than one. Thus 
small absolute changes in the data can be expected to produce large absolute 
changes in the values of u,. 


In Section 3 only the effect of small changes in u, were examined. In the 
mortgage problem the small change was in the value of M. This illustrates 
the fact that a problem which is absolutely ill-conditioned will usually be 
sensitive to small changes in any of its data. 


However, we do not discard this model even though it is absolutely ill- 
conditioned. The over-payment of £427 can easily be sorted out by the 
building society by adjusting the payments in the final year, for example. We 
do have a reasonable model to work with. 
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There were two other examples of how absolute changes in the data gave rise 
to large absolute changes in the solution. 


(1) Increasing the interest rate to 15.25% so that the coefficient of u, 
changes from 1.15 to 1.1525 alters the value of u,; from — 425 to 9833. 
(11) Increasing M from £258 to £263 alters u,, from 6788 to — 1471. 


This last example means that paying £5 a month extra would result in 
approximately 40 fewer payments! 


In the pre-programme exercise you should have calculated that the problem 
of using the recurrence relation to calculate the required values is relatively 
ill-conditioned for values of u, near £20640. Since the starting value used is 
£20000 this is close to the critical value. Hence we would expect small 
relative changes in the data to have a significant effect on the solution. 


If you want to go into this a little further, notice that in the solution of the 
recurrence relation 


U, + ; = au, = p 
given on p. 26 as 


ae P ate’ 
t= [uy + 2 a-—1 


the behaviour of the solution (decreasing or increasing) depends on the sign 


a small change in the data (uo, p or 


E ) so that if uy ~ 
a—1 


a) can change a decreasing solution into an increasing solution. This 
qualitative change in the behaviour of the solution caused by the relative ill- 
conditioning shows up dramatically as we increase the interest rate by small 
amounts. 


of [1 “| 


U, 16% 15.5% 
15.48 
20 000 ak 
ao) 
o 
3 
ie) 
E 
2 We have joined 
Z up the points to 
10 000 Show the shapes 


of the curves 


0 10 20 30 40 Year n 


As we change the rate from 15.4% to 15.5% the solution changes from a 
decreasing solution to an increasing solution. The interest rate of 15.48 % 
corresponds to the interest rate when we simply pay off the interest each year 
but never pay off any of the original loan. If the interest rate increases 
further the repayments are not sufficient even to pay off the interest and we 
end the year owing more than we did at the start of the year. 


Exercise 1 

What would be the monthly repayment for a £20000 mortgage at 15% over 
(i) 15 years 

(11) 20 years 

(111) 30 years 

(iv) 35 years? 

Comment on your results. 


(Solution on p. 55) 
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Exercise 2 
If the repayment on a £20000 mortgage at 15% is 


(i) £263 
(ii) £270 
per month then how long does it take to pay off the mortgage? 


Hint: To find n in the equation a" = b where a and b are known, we take logs of both 


log b 
sides to give nloga = logb or n= = 
loga 


(Solution on p. 55) 


Exercise 3 

If the repayment on a £20000 mortgage is fixed at £258 pounds, how long will it take 
to pay off the mortgage if the interest rate is 

(i) 12% 

(ii) is 

(iii) (14% 

(Solution on p. 55) 


4.2 Coping with ill-conditioning 


In the mortgage problem we saw an example of a problem which was ill- 
conditioned and yet the ill-conditioning could be tolerated because we knew 
how to live with the difficulty. In problems such as the integration problem 
in Section 3, Example 2 the ill-conditioning clearly cannot be lived with. 


There is no method of obtaining satisfactory solutions to this problem. The 
way of overcoming the ill-conditioning here is to change the problem. If we 
can reformulate the problem as a well-conditioned problem we will be able 
to compute satisfactory solutions. 


For the purposes of this discussion we re-examine the problem of finding, for 
example, the integral I,, where 


1 
i= | xe ax 


0) 


Our first approach was to reformulate this problem as a problem involving 
recurrence relations as: 


Find I,, 


using |, = }-—ri,_, Formal statement of the problem. 


with I], =1—e™' 


We saw that this problem is highly ill-conditioned (both relatively and 
absolutely) because a small change in the initial condition ‘wrecked’ the 
calculations. 


Suppose that we were given a starting value other than a value for I). For 
example, suppose we knew that I,, = 0.0455. We could utilize this 
information to work backwards to I,, by turning the recurrence relation 
around to give 
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This is the same recurrence relation as above but we have now expressed 
I,_, in terms of J,. I could now reformulate the problem as 


Find I, 


using), = Reformulation of the problem. 


with I,, = 0.0455 


The method of working backwards using a recurrence relation is known, 
rather unsurprisingly, as backward recurrence. The following table shows my 
calculations 


The first thing to notice about the above figures is how much more plausible 
they are than the results in Section 3. A rather more suprising discovery is 
that if you carry on working backwards to J, you will obtain the computed 
value I, = 0.63212056 which is correct to eight decimal places! This 
reformulation has given us a problem which is very well-conditioned indeed. 
The reason for this can be hinted at if we look at what happens if I, is 
changed to J, = I, + ¢ where ¢ is a small error. If we substitute this into the 


recurrence relation 


1-f, 
eae = r 
we obtain 
a iI, approximation to 
r : 
oe L, , given that i 
= : Z = IS not Known exactly 
ast re4 mae S 


Hence the error in J,_, caused by an error ¢ in J, is 


= E 
I e= 4 I oes Pages r 

Since r is greater than one, this error has been reduced. We can conclude 

that if we begin using I,,, for example, then, using backward recurrence, any 

error in I,, will be rapidly reduced as we compute I,,, I,3, 1, and so on. In 

fact, if the error in J,, is ¢, the error in I,, is —¢/20 and the error in I,, is 

— (—e/20)/19 = ¢/380 and so on. By the time we have computed I,, the 

error ¢ in I,, has been reduced by a factor of 27 907 200 and would be fairly 

insignificant even if ¢ were large. 
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You may have accused me of cheating when I suddenly came out with the 
value 0.0455 for I,, as there was no obvious way of obtaining this value 
easily. However the above remarks indicate that we don’t really need to 
know I,, very accurately in order to obtain accurate approximations for I, 5. 
From our observations of the graphs of the integrands in Figure 1 of Section 
3 we at least know that I, lies between 0 and 1. The following table gives 
the sequence produced by the recurrence relation for each of these values of 
I,). The values for I,, = 0.0455 are repeated for comparison. 


£G,, =0) I, (Ino = 0.0455) 


20 0 0.0455 

19 0.05 0.047725 0 

18 0.05 0.05011974 0.05263158 

17 0.05277778 0.05277113 0.05263158 

16 0.05571895 0.05571935 0.05572755 

15 0.05901757 0.05901754 0.05901703 

14 0.06273216 0.06273216 0.0627322 

13 0.0669477 0.0669477 0.0669477 
0.07177325 0.07177325 0.07177325 


As you can see, whatever starting value we have for I,, between 0 and 1, the 
values for I,, and I,, will be the same correct to eight decimal places. 


Can all ill-conditioned problems be reformulated as well-conditioned 
problems using backward recurrence? The answer is that whenever we have 
a problem of the form 


U4 4 = a,u, + P, 


and |a,| > 1 for all values of r, backward recurrence will help. This is because 
in backward recurrence we have 
1 D, 


r ea a 


r r 


and any error in u,,, will be reduced in the computation of u,. The only 
difficulty in using backward recurrence is in obtaining a suitable starting 
value. : 


Exercise 4 
(i) Show that the following problem is relatively ill-conditioned and reformulate it in 


a form suitable for backward recurrence if you know that uo lies between 0 and 1. 


Find u,, 
using u.,, = 10u, — 3 
1 

ith u, = = 

with uy 3 
(ui) Use backward recurrence to solve for u,. for each of the following initial 
conditions: 

(a) Uy =1 
(b) uy =0 
(C) Uz = 100 


(Solution on p. 56) 

Exercise 5 

The general solution for the recurrence relation 
t.,= 14-5 

5 


a ee 
— a 
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(i) Use your calculator to compute u,;, given that uy = 5/36. 
(ii) Use backward recurrence to compute u,., given that u,, = 21. 
(111) Compare your results (i) and (ii) with the correct values for u, «. 


(Solution on p. 56) 

Exercise 6 

The recurrence relation 
U, + Z — 1 + r 7, 


is to be used to generate a sequence of numbers. If it is known that the terms in the 
sequence tend to zero as r increases, determine an accurate value for u, by using 
backward recurrence with u,,. = 0. 


What is the calculated value of u,, if you start with your calculated value for u, as 
the initial condition? 


(Solution on p. 56) 


Summary of Section 4 


1. The recurrence relation representing the mortgage problem, in which £X 
is borrowed at I % interest and where the monthly repayment is £M, is 


I 
i= { +)" 12M 


The initial condition is uy = X. 
Mortgage problems are often highly sensitive to small changes in data. 


2. Ill-conditioning involving the first order recurrence relation 
Uy = a,U, 52 P, 
with |a,| > 1 for all values of r, can be cured by reformulating the problem as 


a problem involving backward recurrence using the recurrence relation 


— Us =, 
_ 
a, 


A suitable starting value must be obtained in order to use this recurrence 
relation. 


5 Further exercises 


Exercise 1 
Classify the following recurrence relations by completing the table: 


constant- | homo- 
order? | linear? | coefficient | geneous? 


(i) U4 =. ru, 

= 1 7 

il 25 7S 

( ) ret - [, 4 
(iii) i 34, ee 
(iv) Vr+1 = 3ry, + Ty, —2 
(v) Ze+1 a z + 22,2, —1 


(Solution on p. 57) 
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Exercise 2 
Suppose that in a protected environment, in which there are no deaths (!), female 
seals reproduce according to the following rules. 


(i) Each productive female alive at the start of a year produces, on average, 
0.75 new female seals in that year. 


(11) A new born female seal becomes productive in her second year and is 
capable of producing new seals in each subsequent year. 


Express this problem as a problem involving a recurrence relation, given that there 
are 100 female seals at the beginning of year 0 and 120 at the beginning of year 1. 


(Solution on p. 57) 

Exercise 3 

Find the general solution for the recurrence relations 

(i) U,+4 = 7u, se 2 

(ii) u., =u, + 3. 

Find the particular solution of each recurrence relation with u, = 3. 
(Solution on p. 57) 

Exercise 4 


Find the general solution of the linear homogeneous second order recurrence 
relations 

(i) a, =u + O796,_, 

(ii) i, = Sa, = 10m 


Find the particular solutions for each of the recurrence relations for the initial 
conditions u, = 100, u, = 120. 
(Solution on p. 57) 
Exercise 5 
Use your calculator to compute u,, u,, U3, u, and u, using the recurrence relation 
u 
ti, = — +1 
eee 
with (i) uy=l 
and (ii) wu, = 1.0001. 


Comment on your solution. 


(Solution on p. 57) 


Exercise 6 

Use your calculator to compute u,, u, u3, u, and u, using the recurrence relation 
Ung, =(r+1)?u, +1 

with (i) uy = —1.28 
(ii) ut = —1279. 

Comment on your solution. 


(Solution on p. 58) 


Exercise 7 
Given that all the terms in the sequence u;, Uz, u3,..., generated by the recurrence 
relation in Exercise 6, lie between 0 and —1, find u, correct to five decimal places. 


(Solution on p. 58) 


6 Computing 


6.1 Computer packages 


To help you to understand some of the numerical ideas and difficulties which 
arise from using the computer to solve mathematical problems we have 

developed computer packages which are to be used in conjunction with some 
of the units in this course. It is intended that, as part of your work for Units 
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1 and 2, you should make one visit to your nearest computer terminal to use 
the computer to explore problems which would be tedious and very time- 
consuming if tackled using a calculator. 


The computer packages are computer programs, written by Student 
Computing Service programmers, which are designed so that you can tackle 
specific mathematical tasks. The detailed specifications of the individual 
packages can be found in the Student Computing Service Student’s Guide 
check the supplementary material for notification of any modifications. None 
of the packages require any programming skills as they allow you to input 
and solve your problems by answering a series of questions posed by the 
package. The package will also attempt to help you if you get stuck. 


The detailed instructions for logging on to the computer and using the 
packages can be found in the Student Computing Service Student’s Guide 
which has been included with this mailing. If you are not familiar with the ~ 
material in this guide you should read it before you make your terminal visit. 
Note that the packages are held in the program library for this course and 
you can obtain them by following the instructions for ‘Using a Library 
Program’ in the Guide. 


When you have logged on, obtained the program and started to run it, the 
message 


OPTION? 


will be typed at your terminal. This indicates that you have to choose what 
to do next. Your response to the question OPTION? will usually be to 

(i) set up the problem or change it, 

(ii) select the method of solution, 

(iii) solve the problem or stop the program, 

(iv) request help if you are stuck. 


After carrying out any task the program will return to this option point. 


To show you how to use the program, we include an example of how to 
input a recurrence relation using the computer package RECREL. The 
detailed descriptions of the options will be found in Unit 2. 


Example 1 
~ We wish to look at the first ten terms in the sequence generated by the 
recurrence relation 


“44 = = (e+ 1)u, 


with u,. = 0.63212056. (This is the recurrence relation used in Example 2 of 
Section 3.) The following terminal dialogue shows how I input this problem, 
after logging in and obtaining the library program RECREL. Each of my 
responses is underlined to indicate that this is what I typed. The other 
characters in the dialogue were output by the computer. 


44 MST204 1 6.1 


Terminal Dialogue 


OPTION? 10 
ORDER OF THE RECURRENCE RELATION (TYPE 1 OR 2) 
ORDER =?1 


OPTION? 11 
INPUT THE RECURRENCE RELATION 
U(R +1) =?1-—(R + 1)*U(R) 


OPTION? 30 
INPUT INITIAL CONDITION 
U(0) =? 0.63212056 


OPTION? 35 
INPUT NUMBER OF TERMS 
N=?10 


OPTION? 20 
FORWARD RECURRENCE 


OPTION? 41 
FULL PRINTOUT 


OPTION? SOLVE 


The specification of each of the options is contained in Unit 2. 


Comments 


Option 10 enables you to tell the computer the 
order of the recurrence relation. 
The recurrence relation is of first order. 


Option 11 enables you to tell the computer your 
recurrence relation. 

This stands for u,,, = 1 — (r+ 1)u, and should be 
input in this ‘standard form’ (see Subsection 6.3). 


Option 30 enables you to tell the computer the 
initial condition. 
The first term of the sequence must always be up. 


Option 35 enables you to tell the computer how 
many terms in the sequence you wish to output. 


The first ten terms in the sequence are required. 


Option 20 tells the computer the method you wish 
to use (forward recurrence). 


Option 41 specifies what sort of printout you 
require (full printout). 


Option SOLVE tells the computer to perform the 
calculations. 


The 


advantage of this method of input is that you can change any part of the 
problem without having to start again. For example, if you wish to look at 
the first fifteen terms in the sequence you would only have to change N using 


Option 35. 


Notice that for this example I needed the following information: 


(i) the order of the recurrence relation 

(ii) the recurrence relation itself 

(iii) the initial condition 

(iv) the required number of terms in the sequence 
(v) the method to be used 

(vi) the type of printout required. 


By looking at the list of options in Unit 2 I discovered that the appropriate 
options were 10, 11, 30, 35, 20 and 41 respectively. It does not matter in 
which order I give this information to the program but it all has to be given 
before I use the option SOLVE. If you forget to input any of the information 
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when you try to SOLVE your problem, the program will tell you what 
information is missing and then return to the option point where you can 
input the missing data. 


For each package there are five categories of option which can be used in 
response to the question OPTION? 
(1) Command Options 


The command options are specified by name. The following commands will 
be found in each of the packages 


OPTIONS —tthis tells the program to print out a list of the available 
options for the package 


SOLVE  -—this tells the program to run the problem 

HELP —this tells the program that you are stuck and need 
advice. 

LIST —this tells the program to list the information on the 


problem you have input. 
STOP —this is the only way to stop the program. 


Other commands may be used in some packages. 


(ii) Problem Options 


Options numbered between 10 and 19 are used to submit problems or to 
change problems. 


(iii) Method Options 


Options numbered between 20 and 29 are used to choose the method used 
to solve the problem. 


(iv) Data Options 


Options numbered between 30 and 39 are used to submit data such as the 
number of terms, initial values and so on. 


(v) Print Options 


Options numbered between 40 and 49 are used to decide how much 
printout you require. 


HELP 


As stated earlier the command option HELP can be used to obtain guidance 
on how to proceed. However, if you want specific help with a particular 
option, which requires data to be input, you can ask for assistance within the 
option. Here is an example of the HELP command being used within an 
option. 


Example 2 
Suppose that you are not sure how to input a recurrence relation when using 
the package RECREL. To obtain help you would proceed as follows: 


OPTION ? 11 
INPUT THE RECURRENCE RELATION 

U(R + 1) =? HELP 

YOU SHOULD INPUT AN EXPRESSION IN STANDARD FORM TO 
REPRESENT THE RIGHT HAND SIDE OF THE RECURRENCE 
RELATION USING U(R) TO DENOTE THE PREVIOUS TERM AND, IF 
NECESSARY, U(R — 1) TO DENOTE THE PREVIOUS TERM BUT ONE. 
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AN EXPLANATION OF STANDARD FORMS CAN BE FOUND IN 
SECTION 6.3 OF UNIT 1. THE FOLLOWING ARE EXAMPLES OF VALID 
INPUTS. 


U(R + 1) =? 3«U(R) — 4+#U(R — 1). 

U(R + 1) =? 3«SIN(R)*U(R) + U(R — 1) ft 2 
INPUT YOUR RECURRENCE RELATION NOW 
U(R +1) =? 


6.2 The terminal visit 


Before you go 
Before you visit the terminal you should thoroughly prepare what you are 
going to do. You should determine 


(1) how to log on to the computer and run a library program. (This 
information can be found in the Student Computing Service Student’s 
Guide.) 


(11) what problems you wish to tackle. (A list of problems for Unit 1 can be 
found in subsection 6.4.) 


(111) how to input and solve each problem. (Check which options you need 
for each problem after reading Section 6 of Unit 2.) 


If you have a clear idea of what you are going to do you will make optimum 
use of your time at the terminal. 


The visit 

When you make your visit to the terminal you should take with you 

a copy of the Student Computing Service Student’s Guide plus the - 

relevant units. For the package RECREL you will require Units 1 and 2. For 
each of the computing units there is a separate sheet which lists the available 
options. Take this with you and use it to find the correct option if you need 
to choose options while sitting at the terminal. 


The postal service 

If you find it impossible to visit a terminal a postal computing service is 
provided. This is not recommended except in extreme cases since the 
advantage of using the package yourself is that you can react to the output 
from the computer, changing the data, the method of solution and so on. 
Information on how to use this postal service can be found in the Student 
Computing Service Student’s Guide 


6.3 Valid input expressions 


In several of the packages you will be asked to input an expression which 
contains variables, functions and arithmetic operators. There are rules for the 
input of expressions which ensure that the computer interprets the expression 
in an unambiguous way. These rules are almost identical to the rules 
governing the evaluation of an expression such as x”? + 3x + 9, given a value 
of x, using the calculator. 


One way to input the expression x* + 3x + 9 is to type in 


% stands for ‘multiplied by 


XT2+34#X+9 4 stands for ‘raised te the power 


and this would be interpreted correctly by the computer for any value of X. 
This example illustrates some of the rules for arithmetic operators. Note that 
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the computer will evaluate 3 « X before adding it to XT 2. We say that 
multiplication has a higher precedence than addition. Similarly we would 
evaluate Xf 2 first before doing any additions. 

If we wished to input x ‘7*») we would have to use brackets and type 
X7(24+ Y) 


to clarify the meaning. 

Anything in brackets will be evaluated first so that 
(X—Y)T3+(X+Y)T4 

stands for (x — y)> + (x + y)*. 


The order of precedence for the five arithmetic operators is as follows 


level 1 ( ) (brackets) 


level 2 1 (raised to the power) 
level 3 * (multiplied by), / (divided by) 
level 4 + (plus), — (minus) 


For example the computer evaluates an expression such as 
(X+Y)T3*Z4+2*«X/Y 
in several stages: 


(i) evaluate the expression inside the brackets first ie. (X + Y) 

(ii) carry out the operations using f ie. (X + Y)T3 

(iii) carry out all operations using * and /ie. (X + Y)T3* Zand 2+X/Y 
(iv) carry out all operations using + and — ie. (X+ Y)T3*Z+2+#X/Y 


Hence the expression represents 
ox 
(x+y)Pz+ = 


There is a final rule needed to resolve the problems such as 
2/Y¥ *X 

and 

oY 4 


The first could be interpreted either as 2/yx or 2x/y depending on the order 
of the operations carried out. Similarly the second expression could be 
interpreted either as (2 — y) + x or 2 — (y +x). To avoid this ambiguity of 
operators at the same level of precedence we make the rule that operations 
at the same level are carried out from left to right. 


z 2 
Hence 2/Y « X represents : xx = 7 


and 2 — Y + X represents (2 — y) + x 


Exercise 1 
Interpret the following input expressions as algebraic expressions 
(1) 2T(3*X+ Y)—X*Z/Y 


(11) 2tX#*34+4*X*#Y —(34+2Z)T2 
(Solution on p. 58) 
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Exercise 2 
Express the following algebraic expressions as valid input expressions 


(i) (x + y)? + 3xy — z 
@ yrs -S4+— 


(111) 
(Solution on p. 58) 


Functions 
There are several built-in functions which you can use when you input 
expressions. They consist of a three letter name followed by an expression in 
brackets. The following functions are allowed: 
ABS _ the absolute value of the expression ABS(X) means |x| 
e.g. ABS(3) = 3, ABS(— 3) = 3 
EXP the exponential function EXP(X) means e”*, exp(x) 
eg EXP) =< 
INT the largest integer less than or equal to the value of the 
expression 
eg. INT(4-3) = 4, INT(—4-3) = —5 
LOG _ the natural logarithm (to base e) 
LOG(X) means log,(x) 
SOR _ the positive square root of the expression 
e.g. SOR(4) = 2 
SIN the sine of the expression (in radians) 
COS _ the cosine of the expression (in radians) 
TAN _ the tangent of the expression (in radians) 
ATN _ the arc tangent of the expression. 
7 


ATN(X) is in radians and will lie between 5 and = 


Here is an example of the functions being used 
SOR(X 72 + YF 2) — SIN(X + LOG(Y)) — TAN(X) 


The corresponding algebraic expression is 


./x? + y? — sin(x + log,(y)) — tan(x) 


6.4 Computer exercises 


Remember to plan your work carefully, as described in Subsection 6.2, before 
running these exercises at a terminal. 


Exercise 3 
Use forward and backward recurrence on the recurrence relation 
U,+4 =1—(r+1)u, 


r 


with ty = 1 — eo = 65212056... 


which was used to find the values of the integrals discussed in Sections 3 and 4. 
Compute the first 25 terms. Choose a suitable value for u,, for backward recurrence. 


MST204 1 Solutions 49 


Exercise 4 
The recurrence relation 


“44> (P+ 1)u, ee, 
with uy =e —1 = 1.7182818... 


can be used to evaluate the integrals 


1 

w= | ae ax 
@) 

Use the package to compute the first 25 terms in this sequence. 


Exercise 5 

Check that your mortgage repayments are correct so that you repay the loan in the 
specified time period. (Note that your mortgage may include property and/or life 
insurance which are extras and should be omitted from the calculations.) 


If you do not have a mortgage, work out the monthly payment for a £30000 
mortgage borrowed over 25 years if the interest rate is 17%. Use the package to 
determine the amount outstanding at the end of each year. 


Exercise 6 
Use the package to determine the behaviour of the sequences generated by the 
following second order recurrence relations. 


(1) Ups, =U, —Up-1, Up = 0, u, = 1 

(i1) ieee = 28, Uo = Ou = 1 

(iii) U,+, = 0.9u, — 0.2u,_, + 1500, up = 3600, u, = 4000 
(iv) i, = ee — ton, +4 uw =5 8, = 5 


Appendix: 7 
Solutions to the exercises 


Solutions to the exercises in Section 1 (ii) us, =4, u, = 7, us = 11, ug = 18, u, = 29. 
1. (i) Since u,,, is described in terms of u, and u,_,, it 3. Let u, be the amount of money left in the account at 
is a recurrence relation. the end of r years. In the following year the account 


increases by 10% and is subsequently reduced by the 


(ii) Only u, is described in terms of u, and we have no : 
withdrawal of £4000. We express this mathematically as 


idea what u, is. This is not a recurrence relation. 


(iii) This is a recurrence relation. 


G44 =u, + — 4000 


cS a 
(iv) This is a recurrence relation because it can be used 100 
to generate a sequence. Note however that it is not 

expressed as a formula for u,,,. If we wanted to we could = 1.1u, — 4000 


add one to the subscript r, whenever it appears, to get ne 
The fact that we initially have £20000 in the account is 


| given as uy = 20000. 


hay Oo 
res) Using the recurrence relation we obtain the following 
and this would bring it into the usual form. table 


(v) This recurrence relation requires some manipulation 
to get it into a recognizable form. Reordering the 


equation gives r u, 
2u, +2 =e 3u, +1 ae U,. 
Dividing by two and reducing the subscripts by one gives : poe 
a 2 15800 
eo 3 13380 
4 10718 
5 7789.80 
2. The calculator gives the next five terms as follows. 6 4568.78 
(You may get slightly different results if you use a different 7 1025.66 
calculator.) 
(i). x, = 25, x, = 2.05, x, = 2.006098, x, = 20000001, 
x, = 2. (M101 and MS283 students will recognize, from Since ug is negative there is not enough money in the 
Block 1, Unit 1, that the formula here is that used to find account to be able to withdraw £4000 at the end of the 


the square root.) eighth year. Hence the money runs out after eight years. 
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constant- homo- 
order? linear? coefficient? geneous? 


(i) 
(ii) 
(iii) 
(iv) 
(v) 
(vi) 


(vil) 


(N.A. = ‘not applicable’, for these ideas are not defined for 


non-linear recurrence relations.) 


5. The first five terms in the sequence are 


and we can immediately write down the general solution 
as 


u.=Ax a’ 


6. By putting a = 5 in the formula derived in Exercise 5 
we find that the general solution of 


U, + ee Su, 
u,=Ax 3" 

If u. = 3 then we have A = 3 giving the particular 

solution 


u,=3x 3. 


7. The particular solution can be written as 

u, = a". 
From this we can make the following observations. 
(i) Ifa> 1 the solution grows exponentially. 


(ii) If0<a<1 the solution tends steadily to zero. 


(iii) If a< —1 the solution oscillates while increasing in 
SIZe. 


(iv) If0>a> —1 the solution oscillates and diminishes 
- to zero. 


(v) The only cases we have not considered are those 
where a takes the values 1, 0 or —1. 


If a= 1 then u,,, =u, and u, = 1 for all values of n. 

If a= 0 then u,,, = 0 and u, = 0 for all positive values of 
Nn. 

Ifa = —1 then u,, = —u, and u,, = 1, u,,,, = —1 for 
all values of n. 


In each of the above three cases the size of u,, does not 
change as n increases. 


The following graphs illustrate the behaviour of the 
sequences generated by different values of a, with uy = 1. 


uy, 
ee 
1.0 
0.5 
a=0 
0 ] 2 3 4 5 n 
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(ii) With a = 1.1 and p = 200 we have 
“= Bx (Li)" — 2000 

(iii) With a = 3 and p = 8 we have 
u, = Bx (3)"-—4 


9. Ifa=1 we have 
U,,, =U, + P 

If we put u) = A we obtain 
u,=Uy+p=A+tp 
u,=u,+p=(At+p)+p=A+t+2p 
u,=—u,+p=(A+2p)4+ p=A+ 3p 
U4 =u, +p=(A+t 3p)+p=A+t 4p 

and the general solution is 


u, = A +np. 


10. Using the initial condition to evaluate B in the 
above general solution gives the particular solution 


@ «, = 1000 =8 x (L =2 
B = 1000 
The particular solution is u, = 1000 x (1.1)” 
(ii) u, = 1300 = L.1B — 2000 
1.1B = 3300 
L€. B = 3000 
The particular solution is 
u, = 3000 x (1.1)” — 2000 


(iii) u, =Bx 37-4=5 
ie. 9B = 9 giving B = 1 

The particular solution is 
u, = 3"—4 


11. There are two cases to consider: 


i ass 
We have See a 
a—il 
Pp 
N = B— 
OW Uy ee 


hence B= uy +—— 


8. (i) Using the form of the general solution given by —1 


P Pp P 
= Ba" — = ——]a" — ; 
u, Jae and u,, lu, +2 a Sse 
with a = 1.1 and p = 0 we have (ii) a=1 
u, = B x (1.1) We have u, = A + np. 
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Now u, =A 
hence A = uy 


and u, = Uo + np. 


Thus we can express the solution u, as a function of uo, a, n 


and p as 


Solutions to the exercises in Section 2 


1. (i) The auxiliary equation is 
x* = 7x — 12 
which has solutions x = 4 and x = 3. 
The general solution is 
u, =A x 4"+ Bx 3” 
(ii) Here we have the special case in which we have equal 
solutions since the auxiliary equation 
x* = 6x —9 
has solutions 3 and 3. 
The general solution in this case is 
u, = (A + Bn)3” 
(iii) The auxiliary equation is 
ext 4 
Using the formula for the solution of a quadratic, we have 


147144 14/5 
-=y = 


x 
Ze 


The general solution is 


uA) +3 


2 


i 


(iv) The auxiliary equation is 
x? = 0.9x — 0.2 
ie. x? —0.9x + 0.2 =0 


The formula method gives 


0.9+,./0.81-08 09+ ,/0.01 
Se ee 
Oe + Ol 
ee Ses 


Hence the two roots are 0.5 and 0.4 giving the general 
solution as 


u, = A(0.5)" + B(0.4)" 


x 


2. (i) The general solution is 
u, = Ax 4" + Bx 3” 

Uy = 2 gives A+ B=2 

u, =7 gives 44 + 3B =7 

This has the solution A = B = 1 

The particular solution is 


u,= 4" + 3” 


(ii) The general solution is 
u, = (A + Bn)3" 

Uy = 2 gives A =2 

u, =7 gives (A+ B)3 =7 


B=- 
or 3 
5299 
ae 3 
(iii) The general solution is 


a= af 3] rape) 


2 
Up = 1 gives A+B=1 


its i=d3 
u, = 1 gives a{ts2| +a] =f 


es 


Multiplying the first equation by and then 


subtracting the second gives 
= 1 
a2) af 55) 
—1 
BY = 


, 8 i. J 


The particular solution is 


w= B+} (A) FA) -) 


mee Gs eee 2 2 

(iv) The general solution is 

u, = A(0.5)" + B(0.4)” 
Uy = 2 gives A+ B=2 
u, =7 gives 0.54 + 0.4B =7 
This has solution B = —60, A = 62 
Hence the particular solution is 

u,, = 62 x (0.5)" — 60 x (0.4)” 


3. The general solution for the recurrence relation 


U,+4 icc U, - u,—4 


is given by 
1+. /5\" 1 —./5\" 
= A B 
1 5 1—/5 
Now —— = 1.618 while = = 0.618 
5 n 
Thus as n increases a | increases exponentially 


1 — ./5\" 
while = oscillates and decreases. We can 


conclude that for large n we have 
1+ ./5 
2 


For example, if we use the result for Exercise 2(iii) to 
compute u, then we have 


ee eis 
5+} : ) = 129846 


u, = 4| ) (n large) 


while 
JS (i= dae 
1 5 
and already the term involving ty is the most 


significant part of the solution. 


4. The auxiliary equation for this recurrence relation is 


3 
2 
ee 1 
xX ere 


Le. 
3 
t - Yafee 
x a 


The roots of this equation are 2 and —4 so that the 
general solution of the recurrence relation is 


If you consider the two terms of the solution separately 
you can see that the term A x 2” is doubling when you 
increase n by 1 while the term B x (—4)" is halving in size 
and oscillating. Thus as n increases the second term dies 
out and the solution is approximately 


u, = Ax<2" {nm large) 


There is one exception to this result and that is the case 
when A = 0. If the initial conditions are uy = 1, u, = —4 
then the particular solution has A = 0 and B = 1 as 


1 


Now as n becomes large the sequence tends to zero. 


5. The auxiliary equation is 
x? = 0.7x + 0.3 


The roots of the auxiliary equation are 1 and —0.3. 
Hence the general solution is 


p, = A+B x (—0.3)" 


Using the initial conditions to determine the arbitrary 
constants A and B we have 


Dy = A+B =100 
p, = A—0.3B = 120 


Hence 
1.3B = —20 
Lc: 
B = —15.3846... = — 15.38 to the nearest penny 
and 


A = 100 — B = 115.3846... = 115.38 to the nearest 
penny. 


This particular solution is 


p, = 115.38 — 15.38 x (—0.3)" 


As n increases we can see that the second term in the 
solution diminishes rapidly and p, tends to £115.38. After 
only six terms in Example 1 we had p, = 115.37. 
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Solutions to the exercises in Section 3 
1. Using the formula 


n= (t+ 2 Ja — 
with a = 10, p= —3 and u, = 4+ «€ we have 
1 3) i=3 
= te i ee eee 
c (fp +4) +] 10-1 
1 


=ex 10"+-= 
EX = 


The change in the value of u,, is e x 10'* which is much 
larger than the change (¢) in uo. 


2. (i) The scale factor is a” where a = 3 and n = 20. 
Hence the scale factor is 37° ~ 3 x 10° 

Thus the problem is absolutely ill-conditioned. 

(ii) The scale factor is (—4)*7° ~ 9 x 107’ 


Since this is very small the problem is absolutely well- 
conditioned. 


3. If we start with 
U., =7Ta.-—1 
Subtracting the original recurrence relation gives 
441 — Ua, = ru, — 4U,) 
=rr—1)(u_,—4u,_,) using u, — u, 


=(r—1)@,_, — 4-1) 


= 1!(Up — Uo) 
Putting n = r+ 1 we have 
U, a (n <7 1)! (Ug eas Uy) 


Hence the scale factor is (n — 1)! and the problem is 
absolutely ill-conditioned for large values of n. Since the 
scale factor does not depend on the initial condition u, 
absolute conditioning does not depend on u,. However, as 
you will see in Subsection 3.3, relative conditioning does 
depend on up. 


4. The exact values for u,, are 


(i) 6666 666 667 
(ii) 6665 666 667 
(111) 6667 666 667 
(You may get slightly different results on your calculator.) 


Clearly the small changes in the data have produced very 
large absolute changes in the solution. For example 
changing the initial condition from u, = 1 to u, = 1.0001 
has given rise to a difference of one million in the value of 
U9. The problem is clearly absolutely ill-conditioned. 
However the fact that the first three figures in the solution 
are correct means that a small error in the initial 
condition does not really affect the size of the solution. 
This is discussed further in Subsection 3.3. 


5. The general solution of the recurrence relation 
is 


1 
ee 
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and the initial condition u, = + implies that A = 0. Since 
the scale factor is (—3)*° errors will grow rapidly and will 
quickly swamp the solution which is u,,. = 5. Hence the 
problem is relatively ill-conditioned. 


To verify this we could put %) = 5 + ¢ in which case 


1 1 
iyo = (—3)°e + 5 = 59049 + 5 


Uo — Ui _ 59049e 


= 118098¢ 
U1 1 
s 2 
—S  S 
Uo 1 
2 


Hence the relative error has been magnified by a factor of 
59049. 


6. With u, = 6 we have u,, = 6 
With u, = 6 + ¢ we have ti,, = 6 + € x (5)'° 


Hence the absolute error in u,, is 


= 1 10 
Fi0 - Wig. oe () 


= 0.00097656¢ 


The relative error in U,, is a 0.00016276¢ 
10 


—— ae 
The absolute error in u, is ¢ and the relative error is -. 
6 


Both the absolute and relative errors have been reduced 
considerably and this agrees with the fact that the 
problem is both absolutely and relatively well-conditioned 
since |a| < 1. 


7. Using the calculator with the given initial conditions 
gives the following sequences 


(You may get slightly different results on your calculator.) 


0 0.166667 
1 ; 1.33335 
< 0.1425 59.6675 

3 0.125 2976.375 

4 148811.75 

> 7440580.5 


(i) The first problem is both relatively and absolutely ill- 
conditioned as the scale factor is (50)° = 3.125 x 10° and 
the true solution is u, = 4. 


(ii) For the second problem the particular solution is 


1 1 
ee es 
iy = 


and u, = 7440476.3 


= 000004 1.4 < 10 
u,  7440476.3 . 


: 0.166667 — : 
o = —___. = §,000002 = 2 x 10°° 


6 


The relative error has increased by a factor of 7 in the 
computation of u, but this could well be an acceptable 
increase in the magnitude of the relative error. This 
agrees with the theoretical result that this problem is 
absolutely ill-conditioned. Although is ‘close to’ 5, we 
may well consider that the problem of using the 
recurrence relation with u, = ¢ is not relatively ill- 
conditioned. This is really in the ‘grey area’ between being 
relatively ill-conditioned and relatively well-conditioned. 


8. The first 11 terms were computed using the two initial 
conditions 


(You may get slightly different results on your calculator.) 


0 — 1.7182818 — 1.7182818 
1 —0.7182818 —0.7182819 
y —0.4365636 —0.4365638 
3 —0.3096908 —0.3096914 
4 —0.2387632 —0.2387656 
5 — 0.193816 — 0.193828 
6 — 0.162896 — 0.162968 
7 — 0.140272 — 0.140776 
8 — 0.122176 — 0.126208 
9 — 0.099584 —Q.135872 
10 0.00416 —0.35872 
11 1.04576 — 2.94592 


This is a very interesting set of figures as you can see how 
the errors are building up. The last three terms in each of 
the sequences are totally different indicating that the 
problem of computing u,, with either of these two initial 
conditions is both relatively and absolutely ill- 
conditioned. We might also conclude that there is a value 
of u, in between —1.7182819 and —1.7182818 for which 
the solution remains small since the first sequence is 
beginning to increase rapidly while the second solution is 
beginning to decrease rapidly. 


In fact the recurrence relation can be used to compute 
values for the integral 


1 
u= -| x= Made 


0 


for which the initial condition is 
Up = 1 —e = —1.7182818... 


Solutions to the exercises in Section 4 


1. The general solution for the 15°% mortgage where the 
repayment is £M per month is 

u, = B(1.15)" + 80M 
Since the initial condition is u, = 20000 the first equation 
in B and M will always be 

B+ 80M = 20000 


The second equation depends on the length of the 
mortgage. 

(i) For a 15 year mortgage we have u,, = 0 

ie. B(1.15)'? + 80M =0 


Solving the two simultaneous equations for B and M 
gives 


20.000 
Ba ee... 
i-ais) 
00 — B 
M = - 285.028 ... 


Hence the monthly repayment is £285.03 
(ii) For a 20 year mortgage, we have 
B(1.15)?° + 80M = 0 
_ ae 
1 — (1.15)?° 
M = £266.269... 
Hence the monthly repayment is £266.27. 


1 50E SZ... 


(111) For a 30 year mortgage we have 
a 
1 — (1.15)°° 
M = 2335).— 
Hence the monthly repayment is £253.83. 


== 306-69... 


(iv) For a 35 year mortgage we have 


2000 
ee ee ee 
1 — (1.15)>5 
M = 251.891... 


Hence the monthly repayment is £251.89. 


The following table summarizes our results. 


length 
of 
mortgage 


monthly 
repayment 


£285.03 


20 £266.27 
Pa £257.83 
30 £253.83 


£251.89 


It can be seen that there is a significant difference in the 
repayment between a 15 year and a 20 year mortgage but 
the difference approximately halves for each subsequent 
increase of 5 years. Notice how small the change is 
between a 25 year and a 35 year mortgage. 


2 (i) If the repayment is £263 the recurrence relation for 
the 15% mortgage is 


u.,, = 1.15u,— 12 x 263 
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Le. 
u,,, = 1.15u, — 3156. 


(One way of determining the repayment period is to use 
the recurrence relation to generate the figures. We present 
a different approach here.) The general solution is 


3156 
1.15 —1 
= B(1.15)" + 21040 
Since the initial condition is u, = 20000 


we have B + 21040 = 20000 
ie. B = —1040 

Hence 
u, = —1040 x (1.15)" + 21040 


We wish to know the value of n when u, = 0. 


u, = B(1.15)" + 


Le. 
1040 (1.15)" = 21 040 
21 040 
: take logs of 
an 
bot i 
| <a oth sides 
oe 
1040 
= ——_—__—_ = 71.5166... 
Tne (1.15) 


Hence the mortgage is repaid after approximately 215~ 
years. This agrees with the results obtained in the 
television programme. 


(ii) _F or a repayment of £270 we have 
u,,, = 1.15u, — 3240. 

The general solution is 
u, = B(1.15)" + 21 600 

and the initial condition u, = 20000 gives 
B = —1600. 

To find the repayment period n we have 
1600 x (1.15)” = 21 600 


_ log (13.5) 
~ Jog(1.15) 


Hence the repayment period is approximately 183 years. 


= 18.622... 


3. One way of tackling this question is to use your 
calculator together with the recurrence relation 


100) ° 


with uo = 20000 in order to obtain the appropriate values 
of n. 


I 
i ( +a — 3096 


We are going to give an alternative method which leads 
to the same results. 


(i) The recurrence relation for a 12% mortgage is 
44 = 1.12u, — 3096 
The general solution is 


3096 
spain.” 
a as 


= B(1.12)" + 25 800. 
The initial condition u, = 20000 gives B = — 5800. To get 
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the repayment period we need to determine the value of n 
for which u, = 0 Le. 


5800(1.12)" = 25 800 

25 800 

5800 

Taking logs of both sides gives 


and (1.12)" = 


Hence the repayment period is approximately 13 years. 
(ii) The recurrence relation is 
u,, = 1134, — 3096 
The general solution is 
u, = B(1.13)" + 23815.38 
and the initial condition gives B = — 3815.38. 
In this case we have 
= 23815.38 
i= STE = 14,9838. 
The repayment period is approximately 15 years. 
(iii) For a 14% mortgage we have 
u,,, = 1.14u, — 3096. 
The general solution is 
u, = B(1.14)" + 22114.286 
and the initial condition gives B = —2114.286 
| prea 
2114.286 


a ce [raza sz 17916... 
. log 1.14 


Hence the repayment period is approximately 18 years. 


4. (i) The problem is relatively ill-conditioned as the 
coefficient of u, is 10 and the initial condition uy is equal 
to the critical initial condition. 


To use backward recurrence we have to guess the value of 
u, for sufficiently large N. A reasonable estimate would, 
for example, be u,,. = 0 since we know that any error in 
Uy, Will be reduced by a factor of 10 at each iteration. 
Formally we have 


Find u,, 


Ur+4 - 3 


using u, = 
g r 10 


with uz90 = 0. 


(ii) The following table outlines the computations for 


each value of uo. 


u,(Ur. = 100) 


100 


19 | 04 0.3 10.3 

18 | 0.34 0.33 133 

17 | 0.334 0.333 0.433 

16 | 0.3334 0.3333 0.3433 

15 | 0.33334 0.33333 0.33433 

14 | 0.333334 0.333333 0.333433 

13 | 0.3333334 0.3333333 0.3333433 
0.33333334 0.33333333 0.33333433 


You can see that even with u,,. = 100 we still have u,, 
correct to five decimal places. 


0.13888899 
0.97222222 
1.8055556 
2.6388889 
3.4722221 
4.305555 
5.1388852 
5.9721961 
6.8053725 
7.6376078 
8.4632548 
92427833 
9.6994829 
7.88963806 
— 9,7253356 
— 138.07735 


OOrANNMNAHWNK © 


ea ek oe 


(ii) For backward recurrence we use u, = 5 


21 


24 20.142857 
23 19.306122 
pe 18.472303 
a1 17.6389 
20 16.805557 
|e 15972222 
18 15.138889 
17 14.305556 
16 13.472222 


12.638889 


(iii) The particular solution we are looking for has B = 0 
giving 


= ook 
2 
75 $ 2a 
— — + — = — = 12,638889 
Hence u, ; 6 + ~~ 3 


Backward recurrence, even with a small error in u,,, gives 
eight figures correct while the original formulation does 
not even give a solution with the correct sign. 


6. We formulate the problem as 
Find u, 
Ur+4 33 1 


a 


using u, = 
r 


with u,, = 0 


0 
—0.01234568 
—0.0158179 
—0.02073098 
—0.02835364 
—0.04113415 
—0.06507088 
—0.11834121 
—0.2795853 
— 1 2I9I6I3 


i) 


meNmwWwW Bh NWN ~) COO = 


On my calculator I computed u, as 


uy = —1.2795853(023) not displayeu 


Using this value for u, we have 


— 1.2795853 
—0,2795853 
—0.11834121 
—0.06507088 
—0.04113412 
—0.02835312 
—0.02071232 
—0.01490368 
—0.04616448 
4.7393229 


—= © COAT A nA fh WN 


As you can see the errors build up slowly at first, but by 
the time we get to u, they are swamping the solution. The 
value for u,, computed using backward recurrence was 
correct to 10 significant figures. 


Solutions to the exercises in Section 5 


constant- homo- 
coefficient? geneous 


order? linear? 


2. We assume that no seals leave the environment either 
by dying or migrating. The basic equation is 


number 
of female seals at | = 
start of year r+ 1 


number number of 
ot female seals 7 + en born female als 


start of year r in year r 
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Let S, be the number of female seals at the beginning of 
year r. Then we have 


number of 
new born female 


5,75 
r+1 =r F | seals in year r 


The number of seals which are productive in year r is 
S,—, and each seal has, on average, 0.75 female offspring. 
Hence we have 


5,4 = 38, + 0755... 
where S, = 100 and S, = 120. 


3. The general solutions of the recurrence relations are 
1 
(i) u,»=Bxi+ 3 


(ii) u,= A+ 3n 
For u, = 3 the particular solutions are 
8 1 
ee res 
@) 5, 3 x7 + 3 
(nu) u, = 3 + 3n 


4. (i) The auxiliary equation is x? = x + 0.75 and the 
two solutions are x = 3 and x = —4. Hence the general 
solution is 


ff af 


The particular solution is obtained by solving 
uy = A+ B= 100 
oF eee : 
=—-—-—=12 
u, pro 0 
giving A = 85 and B = 15. 


Hence the particular solution is 


85 i +15 a 

= = x |—=]. 

u, x [5 ; 

(ii) The auxiliary equation is x” = 8x — 16 and the two 


solutions are x = 4 and x = 4. Hence the general solution 
is 


u, = (A + Bn)4". 
The particular solution is obtained by solving 
uo = A = 100 
u, = (A+ B4=120 
giving A = 100 and B = — 70. 
Hence the particular solution is 
u, = (100 — 70n)4". 


5. (You may get slightly different answers on your 
calculator.) 


1 1.0001 

2 2.0001 

3 2.00005 

“ 1.6666667 1.6666833 
5 1.4166667 1.4166708 
6 1.2833333 1.2833342 
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A small change in the value of uy has generated a 
different sequence. But note that the difference in the 
values for u; is 0.0000009 which is much less than the 
difference in the values for u, (=0.0001). This indicates 
that the problem is absolutely well-conditioned for small 
changes in the initial condition. 


: = 0.0001 Bas caused a 


The relative change in uy of 


0.0000009 
1.2833333 


much smaller. This indicates that the problem is relatively 
well-conditioned for small relative changes in the initial 
condition. Both of these conclusions are consistent with 
the theory which indicates that if |a,| < 1 for all values of 
r then the recurrence relation 


relative change in u, of = 0.0000007 which is 


Un+4 See a,u, + P, 


will be both relatively and absolutely well-conditioned 
with respect to small changes in up. 


6. (You may get slightly different results on your 
calculator.) 


From the tremendous change in the value of u; caused by 
a small change in u, we can see that the problem is both 

relatively and absolutely ill-conditioned for small changes 
in uy. This is consistent with the theory. 


7. Using backward recurrence with the recurrence 


eae. 
relation u, = nee and starting with 
(i) Uy =9 


(ii) u,=—-—1 


we obtain the following results: 


(You may get slightly different results on your 
calculator.) 


1 0 ~~ 

z —0.01 — 0.02 

3 —0.01246914 —0.01259259 
4 —0.01581983 —0.01582176 
7 —0.02073102 —0.02073106 
6 —0.02835364 — 0.02835364 


The two values of u, are identical to 8 decimal places and 
we could be extremely confident in quoting u, as 


u, = —0.02835 (correct to 5 decimal places). 


Solutions to the exercises in Section 6 


1. (i) 20849 _= 
y 


(ii) 2*x3+4xy—-( +2) 


2 @. M+ Y)13 4 32X4¥ —- YTXR 
(ii) Xt9ef¥ +3)12 — Xe VZ4Z7X 
(iii) (X+3)/(Y + 7) + (6 —X)/(YT2) 


There are variations of these answers which are also 
correct. For example, you may have Z/(X* Y) or Z/X/Y as 
the last term in (ii) and you may have omitted the 
brackets round Y 7 2 at the end of (111). 


3, 4, 5, 6 

We are not supplying answers to the computer exercises 
here. You will be able to obtain solutions at the terminal, 
as explained in Unit 2. 
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